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Abstract 

We study the case of a bright (L>L*) barred spiral galaxy from the rich cluster A 3558 in the Shapley 
supercluster core (0=0.05) undergoing ram-pressure stripping. Integral-field spectroscopy with WiFeS 
at the 2.3m ANU, complemented by imaging in ultra violet (GALEX), B and R (ESO 2.2m WFI), Ha 
(Magellan), K (UKIRT), 24/un and 70/im (Spitzer), allows us to reveal the impact of ram pressure on 
the interstellar medium. With these data we study in detail the kinematics and the physical conditions 
of the ionized gas and the properties of the stellar populations. We observe one-sided extraplanar ionized 
gas along the full extent of the galaxy disc, extending ^13kpc in projection from it. Narrow-band Ha 
imaging resolves this outflow into a complex of knots and filaments, similar to those seen in other cluster 
galaxies undergoing ram-pressure stripping. The gas velocity field is complex with the extraplanar gas 
showing signature of rotation, while the stellar velocity field is regular and the isf-band image shows 
a symmetric stellar distribution. We use line-ratio diagnostics to ascertain the origin of the observed 
emission. In all parts of the galaxy, we find a significant contribution from shock excitation, as well 
as emission powered by star formation. Shock-ionized gas is associated with the turbulent gas outflow 
and highly attenuated by dust (A„=1.5-2.3mag). All these findings cover the whole phenomenology of 
early-stage ram-pressure stripping. Intense, highly obscured star formation is taking place in the nucleus, 
probably related to the bar, and in a region 12 kpc South- West from the centre. These two regions account 
for half of the total star formation in the galaxy, which overall amounts to 7.2±2.2 M Q yr _1 . In the SW 
region we identify a starburst characterized by a ~ 5x increase in the star-formation rate over the last 
~100Myr, possibly related to the compression of the interstellar gas by the ram pressure. The scenario 
suggested by the observations is supported and refined by ad hoc N-body/hydrodynamical simulations 
which identify a rather narrow temporal range for the onset of ram-pressure stripping around t^60 Myr 
ago, and an angle between the galaxy rotation axis and the intra-cluster medium wind of ~ 45°. The ram 
pressure is therefore acting at an intermediate angle between face-on and edge-on. Taking into account 
that the galaxy is found ^lMpc from the cluster centre in a relatively low-density region, this study 
shows that ram-pressure stripping still acts efficiently on massive galaxies well outside the cluster cores, 
as also recently observed in the Virgo cluster. 

Keywords: galaxies: evolution - galaxies: ISM - galaxies: intergalactic medium - galaxies: clusters 
general - galaxies: stellar content - galaxies: clusters: individual: A3558 
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1 Introduction 

Both the current properties and the past evolution of 
galaxies are strongly dependent on their environment 
(e.g. Blanton et al. 120051) and mass (e.g. Baldry et al. 



are unable to identify which process(es) are behind 
these transformations, except indirectly by statistical 
comparison of trend s with predictions from cosmolog- 



20061 . Haines et al. 120071 ). The environmen t is a de- 



ical s imulations (e.g. lBalogh et all 2000; H aines et al 
20091 ). What is fundamentally required are direct 



term i nant for both the mor phology-density (jDresslerl . 



1980 : Dressie r et all 119971) and the star formation 



density dButcher k Oemlerl . 1 19841 ; iLewis et all . |2002| ; 



Kauffmann et al 



"smoking-gun" observations of galaxies in the process 
of being transformed via interaction with their envi- 
ronment. The physical processes behind this trans- 
formation can then be distinguished by resolving their 



2004J) relations: late-type, blue, star- i m p ac ts on the gas contents, kinematics and star for- 
forming galaxies are predominant in the field, while ma tion of the galaxies. 

early-type, red, passive galaxies are preferentially found R am pressure resulting from the passage of the 



in galaxy clusters. This suggests that blue galaxies 
accreted from the field have been transformed into 
the passive SOs and dEs found in local clusters. A 
key question in galaxy evolution is therefore, how 
does the harsh cluster environment act to transform 
the spiral galaxies accreted from the field into passive 
SOs and dEs? 

Several mechanisms dependent upon the environ- 
ment have been proposed and investigated in detail 
and all of them serve to kinematically disturb spiral 
galaxies and/or deplete their reservoirs of gas, and 
so quench star formation. These physical processes 
include gravitational and tidal interactions amongst 



galaxy through the hot and dense intra-cluster medium 
(ICM) can effectively remove the cold gas supply 

1999h and thus 



(Gunn k Gott. 


1972 




Abadi et al.. 



rapidly terminate ongoing star formation in cluster 
galaxies. This process could explain the lower star 
formation rates (S FRs; e.g. Balogh et a l. 2000l ) and 
the redder colours ( Bamford et al. . 20091 ) seen in clus- 
ter galaxies with respect to the field population. Nev- 
ertheless, nmijaressure stripp ing, as originally pro- 
posed by I Gunn k Gott] ()1972l ) requires, in principle, 
the presence of a dense ICM. Thus, its evolutionary 
effect would be limited to cluster cores where the gas 
discs of massive spirals are rapidly truncated. 
1996| ) , T n j as t decade ram-pressure stripping (RPS) 
between galaxies and the cluster gravit ational field w h P PT1 extensively studied using hydrodynamical 
(Byrd k Valtonen 1990), galaxy mergers tearnes k Her aa u^|ft ions (e .g. R oediger k Hensler 2005; Roedig 



galaxies (jToomre k Toomrd . ll972l ; lMoore et al 



199J), gr oup-cluster co l lision s (|Bekki|, [2001|), ram- ger k Bruggerl2006l: Kronberger at al. l2008al : Kaoferer 



pressure (I Gunn k Gottl . 1 197211 and viscous stripping 
(jNulsenl . 119821) . e vaporation dCowie k Song"alial .ll977l) 
and 'starvation' ( Larson et all Il980l ). 

Since these mechanisms are characterized by dif- 
ferent time-scales and efficiencies which depend, in 
turn, on the properties of both the galaxies and their 
environment, they can differ entiall y affect the g alaxy 
properties (Boselli k Gavazzi 20061 ; Haines et al. 2007 ). 
Over the last ten years, the development of huge 
spectroscopic surveys such as the Sloan Digital Sky 
Survey (SDSS), plus the availability of panoramic 
far ultra-violet - far infra-red (FUV-FIR) data from 
the GALEX and Spitzer space telescopes have al- 
lowed the impact of environment on star formation to 
be quantified in unprecedented detail. However, be- 
cause they only consider the global properties of each 
galaxy (e.g. stellar mass, star-formation rate), they 



et al. l2009bl : Tonne sen k Bryan [2009J; Bekki 120091 : 



Jachym et al. l2009h describing the RPS process as 
function of galaxy parameters (stellar mass, inclina- 
tion, orbit, velocity, mass of the gas halo, structure 
of the interstellar medium, etc.) and the density and 
structure of the ICM, suggesting a number of ob- 
servables which can be used as diagnostics of its oc- 
currence, intensity, an d stage. With a 3-D hydro- 
dynamical simulation, Marcolini et al. ( 20031 ) found 
that RPS may extend to poorer environments for 
low-mass galaxies which, thanks to their lower escape 
velocities, are easier to strip. RPS also directly af- 
fects the density and temperature structure of hot 
halo gas of galaxies and is relevant i n the starva- 



tion mechanism in clusters and groups (ILarson et al 



198d : iBekki et all I2OO2I : iMcCarthv et~ail 2008) 



The characteristic signatures of RPS are the pres- 
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ence of gas outflows, distortion and ultimate trunca- 
tion of the gaseous disc without corresponding distor- 
tion of_the old st ellar component (e.g. Kenney et al 
20041 : 120081 : l201lh . Such features have been ob served 



al. 1201 ll ) and ever l arger samples of galaxie s (e.g 



in numerous clust er spirals (e.g. Vogt et al. 12004 



Chung et al. l200d ). Ram pressure can also compress 



Monreal-Ibe ro et al. l20ld : Crowl & Kenney 120081 : 
Pracy et al. 120091 . 120121 ) have demonstrated the effi- 
ciency of this tool both to investigate the nature of 
the galaxy emissions and to understand the essential 
physics. 

A key aim of the AC CESS projedQ (PI P. Mer- 
luzzi, see Merluzzi et al. [201CJ, hereafter Paper I) is 
to detect the signatures of galaxies caught in the act 
gas may also be present, whose nature a nd evolution of trans formation, using the unique combination of a 
depend also on the properties of the ISM (jRoediger <fe B raiageeBcale IFS survey with the Wide Field Spectro- 
graph (WiFeS) at the Australian National University 
2.3m telescope together with an unprecedented data- 
set (from far-ultraviolet to far-infrared) on the galax- 
ies in the core of the Shapley supercluster (SSC) at 



and shock the interstellar medium (ISM) enhancing 
the s t ar formation in the inner dis c (jBvrd &: Valtonenl . 
199d : iFuiita fc Naeashinal . Il999h . Tails of stripped 



20081 : iTonnesen & Brvanl . I2OI0I ) . 



Considering both the theoretical predictions and 
the observable effects, it follows that the detection 
of the stripped gas is a key step in order to better 
understand the role of RPS, but also of tidal strip- 
ping, in the transformation of galaxies from their ac- 
tively star-forming phase to their passive phase. The 
stripped gas in cluster galaxies has been detected 
thr ough narrow-ban d Ha imag ing (e.g. Gava zzi et 
al. l200ll : Yagi et al I2007L l20ld : Yoshid a et al. I2OO2I : 
Kenney et al. 120081 . Smith et al. |2I)10]),Jhrough HI 
imaging (e.g. Oosterloo & v an Go rkom 120051 : Chung 
et al . [20071; Haynes et al. 120071 : Koopmann et al . 
20081 ) . in X— ray im ages (e.g. Irwin Sa razin 1996; 



Wang et al . 12004 Sun & Vikhlinin 120051 : Sun et al . 
20061 : l20ld : Machacek et al. 120061 : Kim et al. 120081 ) 
and w ith multi-ban d obs ervations (e.g. Cr owl e t al. 



20051 : Cortese et al. 120071 : Abramson et al. l201ll ). 



The presence of gas outside of the galaxy disc 
could also arise from mechanisms other than RPS, 
such as tidal interactions. Multi-band observations of 
both the tail and the galaxy are needed to distinguish 
amongst the different gaseous components in the tail 
(atomic, molecular, ionized), and high resolution ob- 
servations are required in order to resolve the physics 
of both the gas and the stars in the galaxy. 

Integral field spectroscopy (IFS) allows us to re- 
solve the different spatial components in the galax- 
ies, and so measure the dynamical disturbance of 
the stellar component, discover disturbed gas veloc- 
ity fields, and determine the local enhancement and 
spatial trend of the star formation. Recent IFS obser- 
vations of individ ual ga laxies (e.g. Co rtes e t al. 20061 : 



Crowl & Kenn ev 120061 : Farage et al. l20ld : Jimenez 



Vicente et al. 



20ld : 



=0.048 (IMercurio et al 
HE 



Haines et al 



2011a 



2006; Me rluzzi et all . I2OI0I : 



Smith et all 120071 ). includ- 



ing new narrow-band Ha imaging obtained with the 
Maryland-Magellan Tunable Filter (MMTF) on the 
Magellan-Baade 6.5m telescope. 

Our ongoi ng ACCES S IFS survey with WiFeS 
(Dopita et al. I200I l20ld ) commenced in April 2009, 
just after the WiFeS commissioning. The galaxies 
observed with WiFeS sample different environments, 
from dense cluster cores to the regions where cluster- 
cluster interactions are taking place, to the much 
less populated areas. The whole sample consists of 
24 galaxies assessed as belonging to the Shapley su- 
percluster according their spectroscopic redshift, and 
classified as either star-forming, AGN or composite 



galax ies using AAOmega spectroscopy (jSmith et al 
20071 ). They mostly have intermediate IR colours 



(0.15<f24 Aim /fftr<1.0; a proxy for specific-SFR) be- 
tw een st ar-forming and passive galaxies (Haines et 
a. l2011bl . hereafter Paper III) and K < 13 mag (K < 
K*+1.3). All the galaxies are well-resolved in the op- 
tical images, and often display either disturbed mor- 
phology, such as asymmetry and tails, or evidence of 
star- formation knots. These galaxies have full pho- 



Rich et al. 2010; Sanchez et 



lr rhe project ACCESS ( "A Complete CEnsus of Star for- 
mation and nuclear activity in the Shapley supercluster"), a 
European International Research Staff Exchange Scheme of the 
7th Framework Programme involving the Italian Institute for 
Astronomy and Astrophysics - Astronomical Observatory of 
Capodimonte, the Australian National University, the Univer- 
sity of Birmingham and the University of Durham. 
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tometric coverage which provides complementary in- 
formation on their star formation. In this work we 
present the results for the galaxy - SOS 114732, a 
bright (L^L*) spiral galaxy in the rich cluster A 3558. 

In Sect. [2] we describe the target. The details 
of IFS observations and data reduction are given in 
Sect. [3] and the data analysis is described in Sect. 
The narrow-band Ha imaging obtained using the 
MMTF are presented in Sect. HI We present the re- 
sults in Sect. [6] where the morphology of ionized gas, 
gas and stellar kinematics and the dust attenuation 
map are analysed. The physical properties of the 
gas and the star formation across the galaxy are dis- 
cussed in Sect. [7J All our findings strongly suggest 
the occurrence of RPS as outlined in Sect. [8] where 
the nature of the gas outflow is discussed considering 
also other possible causes as galactic winds and tidal 
interactions. The RPS is confirmed in Sect. [S] by ad 
hoc hydro dynamical simulations. We summarize the 
results and draw conclusions in Sect. [TU1 

Throughout the paper we adopt a cosmology with 
n M =0.S, ^a= 0.7, and H =70kms- 1 Mpc- 1 . Ac- 
cording to this cosmology 1 arcsec corresponds to 
0.941 kpc at z=0.048 and the distance modulus is 
36.66. 

2 The galaxy SOS 114372 

The galaxy SOS 114372, named following the Shap- 
ley Optic al Survey (SOS) ident ification (Mercurio et 



al. 120061 : Haines et al. 120061 ). is a bright (L;>L 



Table [2j The galaxy is amongst the brightest SSC 
members in the infra-red and it is also the bright- 
est SSC galax y in the ultra-violet (UV) (see Table [2] 



Haines et al. l2011al . hereafter Paper II), providing 



evidence of the presence of both obscured and unob- 
scured star formation. The total infra-red luminosity 
of the galaxy is estimated as Lrr R=l.l=fc0 : l-10 11 -^m 



based on fitting the infra-red SED models of lRieke et al. 



(120091 ) to t he Spitzer/MIPS 70um flux and AKARI 



90/xm flux ( Murakami et al. . 20071 ). These data are 



also consi stent with the prio r IRAS lOOj/m flux mea- 
surement ( Allen et al. . 1991 ). Miller ( 20051 ) measured 
a radio flux at 1.4GHz of 9.31 mJy corresponding to a 
luminosity Li.4g_hz=5.3-10 22 WHz _1 , consistent with 
expectations from the FIR-radio correlation in which 
both the FIR and radio emissions come from ongo- 



i ng st ar formation (Paper III) . Following iLerov et al 



(|2008h . we estimated a global SFR=8.53l 2 ;^ M yr" i 
(of which 78% is obscured, Paper II) with the errors 
accounting for the SFR calibration uncertainty. The 
most intense star formation occurs in the heavily ob- 
scured centre, but an intense region of less-obscured 
star formation is also observed in the SW disc (see 
Sect. [63]). 

We estimate a stellar mass of A4*=7.0±0.7-10 10 M Q 
(Paper I). SOS 114372 has dynamical mass within 20 kpc 
radius of A / fd l/ri ~2-10 11 M (see Sect 16. 2 p fr om which 



we can estimate the total halo mass following lConselice et al 



()2005h of X=1.2-1O 12 M . To derive the structural 
properties of the galaxy, we modelled the light distri- 
bution in K band with a Sersic profile f or the bulge 



plus an exponential disc using GALFIT ([Peng et al 



if=11.504±0.015 ) spiral galaxy at redshift z=0.0506 
rtDale et al.l . ll999l ) belonging to the rich cluster A 3558 

(Abell richne s s 4) which has a med ian redshift z=0.0477 the effect ° f strong dust absorption. The disc radial 



20101 ). The choice of the X-band allows to mitigate 



([Smith et all 120041 : IProust et all 120061 s ) and a veloc- 
ity dispersion of lOlOkms -1 ( Proust et al. . 20061 ). 



SOS 114372 is almost edge-on and is located at the 
projected distance of 0.995 Mpc from the centre of 
the cluster A 3558, well within the cluster radius 

r 50 o=l-214±0.044 Mpc (see Sanderson & Ponmanin^f 1 rat i° f is tn ] ed with f spect to the disc 
and in a relatively low-density region (p=0.88 gals 



scale length is r^=6.4 arcsec (~ 6 kpc) which, after 
applying the "d ust co rrection" introduced by Gra- 
ham &: Worley ( 20081 ) becomes r& ~ 5.3 kpc. The 
effective radius of the 'bulge' is r e =1.5 arcsec (~ 
1.4 kpc). The 'bulge' is significantly elongated, with 



arcmin ) as measured from the SOS R < 21 galax- 
ies (Fig. 1, left panel). This density is less than half 
of the peak density observed in the core of A 3558. 
The main properties of the galaxy are listed in 



by 13°, and its Sersic index is n=1.2 indicating that 
it is actually a bar. The disc axis ratio is 0.14, which 
implies an inclination of the disc with respect to the 
line of sight of 82° . The i-T-band bulge-to-disc ratio is 
B/D=0.44, after c orrecting for inclination following 



Driver et al.1 (|2008l ) 
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Figure 1: Left: The position of SOS 114372 is marked (blue dot) in the field of A 3558 (North up, East 
left) where the surface density of R<21.0 galaxies from the SOS is shown in gray scale (from 0.33 to 2.0 
galaxies arcmin -2 ). The r5oo cluster radius is indicated by the red circle. Right: BRK composite image 
of SOS 114372. The two white rectangles denote the WiFeS pointings: #1 centred on the galaxy, #2 
covering the NE side. The two galaxies to the NE are supercluster members: SOS 114493 (2=0. 05) and 
SOS 112393 (z=0.046). SW of SOS 114372 is a faint galaxy, SOS 115228, without redshift measurement. 
The orientation is indicated and adopted for all the following galaxy maps. The white bars indicating the 
orientation correspond to a length of 15kpc. Channels R, G, B are assigned to K, R, B bands, respectively. 
The nuclear region appears to be highly obscured by dust. 
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Table 1: Properties of the galaxy SOS 114372. 



Property 


Value 


Source 


Magnitudes and fluxes 


FUV a (1510 A) 


17.11±0.05 


Paper II 


NUV a (2310 A) 


16.68±0.03 


Paper II 


B h 


15.16±0.04 


Mercurio et al. 2006 


R b 


14.27±0.04 


Mercurio et al. 2006 


K b 


11.504±0.015 


Paper I 


24/mi 


42±2mJy 


Paper II 


70/xm 


0.758±0.038 Jy 


Paper II 


90/zm 


0.837±0.086 Jy 


Murakami et al. 2007 


100/im 


1.2±0.2 Jy 


Allen et al. 1991 


1.4GHz 


9.31 mJy 


Miller 2005 


.fT-band bulge-disc decomposition 


disc radial scale length (rj) 


5.3 arcsec 


this paper 


bulge effective radius 


1.5 arcsec 


this paper 


Sersic index of the bulge 


1.2 


this paper 


bulge-to-disc ratio 


0.44 


this paper 


disc inclination 


82° 


this paper 


Velocities 


rotation velocity at 


200±13kms~ 1 


this paper 


radial velocity c 


826±13kms- 1 


Dale et al. 1999: Proust et al. 2006 


Masses 


Stellar mass 


7.0±0.7-10 iU M o 


Paper I 


Dynamical mass rf 


2-10 n M Q 


this paper 


Total halo mass 


1.2-10 12 M 


this paper 


Distances 


Redshift 


0.0506 


Dale et al. 1999 


Projected distance to cluster centre 


0.995 Mpc 


Sanderson & Ponman 2010 


Star formation rates 


Global SFR from UV+IR 


Moyr" 1 


Paper II 


Global SFR from Ha 


7.2±2.2M yr- 1 


this paper 



a) AB photometric system. 

b) Vega photometric system. 

c) With respect to the cluster systemic velocity. 

d) Within 20kpc radius. 
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The BRK composite image, as der ived from the 
SOS (ESO-WFI, Mercurio et al, l200d ) and if-band 
surveys (UKIRT-WFCAM, Paper I), shows evidence 
of dust absorption in and SE of the galaxy centre 
(Fig. [U right panel). In the composite image we also 
notice hints of matter beyond the stellar disc (north- 
west from the disc, see Sect. [H]). The orientation of 
Fig. Q] with North upwards and East to the left will 
be adopted throughout the article. 

3 Integral-field spectroscopy 
3.1 Observations 

The spectroscopic data of SOS 114372 were obtained 
during two observing runs in April 2010 (pointing 
#1) and April 2011 (pointing #2) usingthe Wide- 
Field Spectrograph (WiFeS; Dopita et al. I2007L lioich 
on the Australian National University 2.3m telescope 
at the Siding Spring Observatory, Australia. WiFeS 
is an image-slicing integral-field spectrograph that 
records optical spectra over a contiguous 25"x38" 
field-of-view. This field is divided into twenty- five 
l"-wide long-slits ('slices') of 38" length. The spectra 
were acquired in 'binned mode', with 25x38 spaxels 
of l"x 1" size. Hereafter we will refer to the spatial 
directions across the slices or along them as 'X' or 'Y' 
direction respectively. WiFeS has two independent 
channels for the blue and the red wavelength ranges. 
We used the B3000 and R3000 gratings, allowing si- 
multaneous observations of the spectral range from 
~3300A to ~9300A with an average resolution of 
R=2900. For further details on t he WiFeS instru- 
ment see Dopita et al. (120071 . l2010h . 

The positions of the WiFeS field-of-view for the 
two pointings are shown in Fig. [T] (right panel). The 
total integration time on the galaxy was of 4.5 h in 
run #1, obtained as the sum of 6x45min exposures. 
Pointing #2 was obtained to map the gas distribution 
in the Northern side of the galaxy and to measure the 
extent of extraplanar gas. The data of this pointing 
are shallower (2x45 min) than run $d and are used 
to map the gas and derive the kinematics, but not to 
measure line flux ratios. For each galaxy spectrum, 
we also acquired the spectrum of a nearby empty sky 
region with 22.5 min exposure to provide accurate sky 



subtraction. 

For each galaxy exposure we obtained spectra of 
spectrophotometric standard stars for flux calibra- 
tion and stars with nearly featureless spectra to mon- 
itor the atmospheric absorption bands. Arc and bias 
frames (see below) were also taken for each science 
exposure. Internal lamp flat fields and sky flats were 
taken twice during both runs. 

3.2 Data reduction 

The data were r educed using the W iFeS data re- 
duction pipeline ( Dopita et al. . 2010l ) and purposely 
written FORTRAN and IDI@ codes. The WiFeS 
pipeline performs all the steps from bias subtraction 
to the production of wavelength- and flux-calibrated 
data-cubes for each of the B and R channels. The 
average r.m.s. scatter around the dispersion relation 
was 0.2-0.3 A. The final spectral resolution achieved 
is cj~40kms~ 1 , and is wavelength and position de- 
pendent (see below). The data-cubes were sampled 
at 1" x 1" x 1 A and cover a useful wavelength range of 
3600-9000 A. The data-cubes are also corrected for 
atmospheric differential refraction. 

Accurate sky subtraction is particularly impor- 
tant for our data because of the presence of OH lines 
in the vicinity of the Ha line and the [Nil] doublet at 
the redshift of the SSC. Sky subtraction was carried 
out as follows. The sky spectrum taken closest in 
time to the galaxy spectrum was first cleaned of cos- 
mic r ay hits using the code LACosmic (jvan Dokkuml . 
200lh to detect them, followed by an interpolation 
along the spatial direction to remove them. The sky 
spectra were then smoothed in the Y direction with 
a median filter. This allowed us to effectively min- 
imize the introduction of noise in sky subtraction. 
In an ideal case, we should have multiplied the sky 
cubes by a factor of two (the ratio of exposure times) 
before subtraction, but due to differences in airmass 
and changes in sky conditions, especially affecting the 
amplitudes of emission lines, we normalized the sky 
frames to the galaxy frames by comparing the relative 
strengths of the atmospheric emission lines. This re- 
sulted in a range of multiplicative factors for the sky 



2 http: / / www.exelisvis.com/language/en- 
US/ProductsServices/IDL.aspx 
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frames (~1.5-2.5). 

Particular attention was devoted to the removal 
of the atmospheric absorption features, since at the 
redshift of the SSC the atmospheric band at ~6870 A 
may affect the [Nil] -Ha group. This O2 band is ac- 
companied by two other O2 bands (~ 6280 A and ~ 
7600 A ) which could be effectively used to monitor 
the quality of the subtraction. 

The co-addition of the individual exposures re- 
quires an accurate spatial registration of the data- 
cubes. For this purpose, we took advantage of the 
high quality ESO-2.2m WFI B- and i?-band images 
of the galaxy from the SOS. The bandpasses of these 
images can be matched to the B and R channels of 
WiFeS. The true spatial scale of WiFeS data was de- 
termined from (i) the distance between the galaxy nu- 
cleus and a star taken during run #1 and (ii) images 
of the centre of the globular cluster NGC 3201 ac- 
quired in run #2 for this purpose. We find that while 
the X scale of WiFeS is l"/pixel, the Y scale is larger 
than the nominal 0.5"/pixel by a factor 1.15 for the 
B data-cube and 1.04 for the R data-cube. After cor- 
recting for these pixel scales we multiplied the WiFeS 
cubes with the WFI response curves and determined 
the relative spatial position by cross-correlation. The 
individual exposures were then registered and co-added 
using IRAF. B and R fluxes were adjusted, where 
necessary, by matching the fluxes in the common 
wavelength range (5300-5600 A ). The adjustment re- 
quired was always within the estimated flux errors 
(see Sect. 15. 2|) . The reduced data-cubes were finally 
corr ected for Galactic extinction following Schlegel et 
al. (19981 ) and using the extinction curve bv lCardelli et 
(|1989l ) with Ry=3.1. The sensitivity of our data 



Chile on 20 May 2012. The MMTF is based on a 
Fabry-Perot etalon, which provides a very narrow 
transmission bandpass (~5-12A) tha t can be tuned 
to an y wavelength over ~5000-9200A (IVeilleux et al 



turns out to be 0.5-10 _17 ergs _1 cm _2 A _1 arcsec -2 at 
a signal-to- noise ratio SNR=5. This is the lowest 
SNR for which we consider our measurements reli- 
able. 

4 Ha imaging 

Deep Ha imaging of the galaxy SOS 114372 was ob- 
tained w ith the Maryland-M agellan Tunable Filter 
CMMTF: IVeilleux et alll2O10h on the Magellan-Baade 
6.5m telescope at the Las Campanas Observatory in 



20ld ). Coupled with the exquisite image quality pro- 
vided by active optics on Magellan and the Inamori- 
Magellan Areal Camera & Spectrograph (IMACS), 
this instrument is ideal for detecting extra-galactic 
Ha-emitting gas. The MMTF 6815-216 order-blockin; 
filter with central wavelength of 6815A and FWHM 
of 216A was used to provide coverage of the Ha emis- 
sion line for galaxies belonging to the Shapley super- 
cluster. 

At the start of the observing run the etalon plates 
were parallelized and the wavelength calibration of 
the etalon for our setup was determined and further 
checked throughout the night. The MMTF does not 
provide a monochromatic image over the full extent of 
the IMACS 27'x27' field of view, but rather a central 
circular monochromatic region (the Jacquinot spot), 
outside of which the central transmission wavelength 
monotonically decreases with increasing angular dis- 
tance from the optical axis in a well understood man- 
ner. The instrumental set-up was chosen to provide 
the largest Jacquinot spot (~11.5') and a transmis- 
sion bandpass of FWHM 10. 3A (corresponding to 
~470kms _1 ). The target galaxy was always placed 
at the same location, well within the Jacquinot spot, 
but away from the gaps between CCDs. 

The galaxy was observed for a total 60 minutes 
in Ha (4x900 sec). The galaxy was also observed for 
15 minutes in the continuum, by shifting the central 
anlavelength of the etalon ~50A bluewards to exclude 
emission from both the Ha line and the adjacent [Nil] 
lines, and into a wavelength region devoid of major 
skylines. The typical image resolution for these ex- 
posures was 0".50. 

These data were fully reduced using the MMTF 
data reduction pipelinqj, which performs bias sub- 
traction, flat fielding, sky-line removal, cosmic-ray re- 
moval, astrometri c calibration and sta cking of multi- 
ple exposures (see Veilleux et al. . 2010l ). Photometric 
calibration was performed by comparing the narrow- 
band fluxes from continuum-dominated sources with 



http:/ /www.astro. umd.edu/~veilleux/mmtf/datared. html 
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their known i?-band magnitudes obtained from our 
existing WFI images. Conditions were photometric 
throughout and the error associated with our ab- 
solute photometric calibration is ~10%. The effec- 
tive bandpass of the Lorenzian profile of the tun- 
able filter of 7r/2xFWHM is then used to convert 
the observed measurements into Ha fluxes in units of 
erg sec -1 cm -2 . The filter is sufficiently narrow that 
there should be little or no contamination from [Nil] 
emission. The data were obtained in dark time result- 
ing in very low sky background levels, with la sur- 
face brightness fluctuations within a 1 arcsec diam- 



eter aperture of 0.2x10 
implying a sensitivity of 



17 



ergs 



_1 cm 2 i 



arcsec 



1.0x10 17 ergs 1 cm 



2 A _1 arcsec" 2 at SNR=5. 



5 Data analysis 

In this section we describe the stellar continuum mod- 
eling and the measurement of the emission line fluxes. 
The emission-line fluxes are used to estimate: i) the 
gas kinematics (Sect. I6.2f) : and in each galaxy region 
ii) the line diagnostics (Sect. \77T\i ; iii) the dust atten- 
uation for the ionized gas and SFR across the galaxy 
(Sects. [6731 and IT72"j) . The stellar continuum modeling 
is needed to obtain the pure emission-line spectrum 
and, accounting for the dust extinction, allow us to 
infer stellar population ages in different galaxy re- 
gions (Sect.ETJ). 

5.1 Stellar continuum modeling and sub- 
traction 

Late-type galaxies present complex star formation 
histories with continuous bursts of star formation from 
the earliest epochs right up until the present day 
dKennicuttlll98.j Kcnnicut t et all 11994 : iJames et al.1 



tended wavelength ranges (i.e. not just absorption 
lines) by linear combinations of multiple stellar pop- 
ulations, while also accounting for the impact of line- 
of-sight stellar motions and instrument resolution, 
and the non-linear effects of dust extinction. A va- 
riety of implementations of full spectral fitting have 
been developed, which explicitly attempt to derive 
robust stellar population parameters and their un- 
cert ainties, including the use of detailed simu l ations 



20081 : IWilliams et all . 1201 ll ). This implies that, un 



like early-type galaxies, their spectra and absorption 
line indices cannot be successfully fit by comparison 
to those from single stellar population models. In- 
stead the technique of full spectral modeling has been 
developed to reliably describe the properties of the 



individual stellar components (|Cid Fcrnande s et al. 



2009). 



Cid Fernandes et all 120051 : iKoleva et aD . 120081 . 



To obtain a reasonable fit to the stellar continuum 
a signal-to- noise leve l of ^> 50/A is generally required 



(e.g. Sanchez et al. 120111 ). We take for each spaxel 



the continuum from the 3x3 region centred on the 
spaxel, which for the main body of the galaxy was suf- 
ficient to achieve the necessary signal-to-noise. Out- 
side the galaxy disc, further smoothing was required, 
and so increasingly large rectangular regions centred 
on the spaxel were considered until a mean signal- 
to- noise level of ~40/A was reached for the stellar 
continuum over the wavelength range 4600-4800 A. 
For each spaxel, the spatially-smo othed spectrum 



from t he blue arm is firstly fit with the lVazdekis et al 



()2O10h stellar population synthesis models. We mask 
out regions below 3950 A which have significantly re- 
duced signal-to-noise levels and flux calibration reli- 
ability, and above 5550 A where a bright sky line is 
located. The wavelength regions affected by emission 
lines are also masked out. 

We consider a total of 40 simple stellar popula- 
tions (SSP) models covering the full range of stel- 
lar ages (0.06-15 Gyr) and three different metallici- 
ties [M/H]=-0.41, 0.0, +0.22. The majority of SSP 
models had solar metallicity and provided the re- 
quired fine sampling of age, while additional evolved 
(1-15 Gyr old) super-solar models were included to 
improve the fitting of the sma ll-scale c o ntinu um fea- 
tures. The models assu me a iKroupal (120011) initial 
mass function (IMF). The lVazdekis et al.1 (|201(ih SEDs, 
based on the Medium resolu tion INT Library of Em 



pineal Spectra (MILES) of ISanchez-Blazquez et al 



(2006), rave a nominal resolution of 2.3 A, close to 
our inst rumental resolution, and cover the spectral 



2005 : lQcvirk et al.l . l2006l : lKoleva et al.l . l2009l : lMacArthuiraetod 540-7410 A 
20091 ) . This involves the fitting of spectra over ex- In a first pass, the best-fitting Vazdekis et al. SSP 
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Figure 2: Example result from the stellar continuum fitting process. The black curve shows the input 
spectrum, coming from the 3x3 spaxels region, centred on the photometric centre of the galaxy. The magenta 
curve shows the resultant best-fit stellar continuum comprising a linear combination of SSPs, requiring both 
old (> 1 Gyr; thin red curve) and young (< 1 Gyr old; thin blue curve) components. The shaded regions 
indicate the wavelength ranges excluded from the fitting process, including the masks for emission lines. 
The residual emission component (thick blue curve) reveals clear emission at [On] A3729, H5, H7, H/3 and 
[Oiii] AA 4959 5007. 
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Figure 3: Example of the stellar continuum modeling in the red arm. The black curve shows the red-arm 
spectrum, coming from the same 3x3 spaxel region as Fig. [2J which is located in the core of the galaxy. 
The magenta curve shows the resultant best-fit stellar continuum comprising a linear combination of SSPs, 
based solely on the fit to the blue arm, and extrapolated to the red arm. The blue curve shows the residual 
emission, i.e. subtracting the magenta curve from the black curve. Notice the emission lines of [Oi], Ha, 
[Nil], [Sn] at the observed wavelengths 6620A, 6895A, 6880-6917A, 7057-7072A. 
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is taken to derive the absorption- line redshift, z a bs, 
and subsequently a first estimate of the velocity dis- 
persion, <7, assumed to be in the range 0-300 kms -1 . 
The value of a can be strongly affected by template 
mismatch, and so it is important that the full range 
of SSP models are fitted over. The a is determined 
as that which minimizes the x 2 value for all combina- 
tions of a and SSP model. We then fix the values of 
z a b s and a, and iteratively fit the observed spectrum 
by a non-negative linear combination of the 40 SSPs 
in order to minimize the x 2 value, until convergence 
is reached. Throughout this iterative process, indi- 
vidual strongly discrepant pixels (>5cr) are identified 
and masked if necessary, to ensure that single pix- 
els do not bias the convergence of the fit. This final 
linear combination is typically found to be composed 
of 8-10 individual stellar populations. A final fit to 
estimate a is made using the best-fit complex stel- 
lar population model. The final linear combination 
of SSPs is then re-scaled linearly to obtain the best- 
fit to the spectrum from the single central spaxel and 
then subtracted to obtain the pure emission spectrum 
for that spaxel corrected for stellar absorption. We 
neglect internal extinction in this continuum fitting 
process, since our aim here is to fit the stellar contin- 
uum in order to model and subtract it in the regions 
of the emission lines, while we include dust extinction 
in our later star formation history analysis limited to 
three galaxy regions (Sect. 17.21) . 

Following Moustakas et al. ( 2O10l ). any remaining 
residuals (typically of the order a few percent), due 
to imperfect sky subtraction, or template mismatch, 
are removed using a 25-pixel sliding median (again 
masking out the emission lines). This last stage had 
little impact on the residual spectrum, except for the 
extreme blue end with A<3800 A, where the flux cali- 
bration accuracy is weakest, and the additional mask- 
ing of emission lines, ensured that it did not impact 
the final emission line measurements. 

Having fit the stellar continuum for the blue arm 
spectrum, this best-fit linear combination of SSPs is 
taken and extended into the red arm. Then vary- 
ing only the global scaling factor to account for any 
(slight) mismatch in the flux calibration between red 
and blue arms, it is subtracted from the red arm spec- 
trum, to produce the pure emission spectrum for the 



corresponding spaxel in the red arm. 

In Fig. [2] we show an example output of our stel- 
lar continuum fitting process. The input spectrum 
(black curve) is from a 3x3 spaxels region in the 
centre of the galaxy. The resultant best-fit stellar 
continuum (magenta curve) comprising a linear com- 
bination of SSPs, requires both young and old (thin 
blue/red curves) components. Almost all of the small- 
scale structures seen in the observed spectrum are 
real rather than noise, being precisely mapped by 
the model stellar continuum. The shaded regions 
correspond to wavelength ranges not considered in 
the fitting process, including the masks for emission 
lines. This stellar continuum is subtracted from the 
observed spectrum to produce the residual emission 
component (thick blue curve) accounting now for stel- 
lar absorption, revealing clear emission at On, H<5, 
H7, H/3 and [Oill]. The Balmer emission lines are 
all located within deep absorption features, demon- 
strating the necessity of accurately modeling and sub- 
tracting the stellar continuum prior to measuring these 
lines. Outside of the emission lines, the rms levels in 
the residual signal are consistent with expectations 
from photon noise, with little remaining structure. 
This holds true throughout the galaxy indicating that 
on a spaxel-by-spaxel basis the model fits to the stel- 
lar continuum are formally good (x 2 ~l)- Figure [3] 
shows the fit extended in the red arm for the same 
spatial region as Fig. The SSP model continuum 
(thick magenta curve) of the spectrum (black curve) 
is able to describe all the absorption features and 
once subtracted (blue curve) allows robust emission 
line measurements. 

5.2 Emission-line measurements 

Emission-line fluxes and widths were measured with a 
purposely written FORTRAN code, which performs 
a Gaussian fit to the emission lines. Where lines 
are either partially overlapping or close (as for the 
groups [Nil] -Ha- [Nil], [Sll] 6717-6731 and [Olll]4959- 
5007), the lines are fitted simultaneously. The fits 
are however left completely independent even in the 
case when their flux ratios are known to be fixed. 
This was done to assess the reliability of our mea- 
surements (see below). The code first fits a straight 
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line to the residual continuum (after sky subtraction 
and stellar population fitting and subtraction) on ei- 
ther side of the line, and then fits a Gaussian to each 
continuum-subtracted line to derive the position of 
the maximum, the peak amplitude, and dispersion 
a. These values are then used to perform a series 
of (50) simulations in which random scatter derived 
from the continuum fit is added to the solution, and 
these simulated lines are fitted again. The standard 
deviation of the distribution of the fitting parameters 
is then adopted as the contribution of the fit to the 
uncertainty. Fluxes are computed from the measured 
amplitudes and a. 

We estimate the SNR of the emission as being the 
ratio of the peak amplitude to the standard deviation 
of the surrounding continuum. Comparing the SNR 
with the error coming from the fit and taking into 
account the appearance of the lines, we fix at SNR=5 
the lower limit for reliable measurements. This limit 
corresponds to ~30% relative errors from the fit for 
both the flux and a. The total error on the flux 
is computed adding in quadrature the uncertainty in 
the flux zero-point, estimated comparing the fluxes of 
the individual exposures taken in good atmospheric 
conditions. This amounts to ~7% -10% (R and B 
respectively). 

Figure H] shows examples of both real and sim- 
ulated data, and the computed fits in the case of 
SNR~5 (our lower limit) and SNR~50 (the average 
of our data) in the wavelength region around Ha. 

The instrumental contribution to the line widths 
was estimated by measuring the widths of selected 
arc-lamp lines from the B and R data-cubes. Due 
to the combination of spherical aberration and high- 
order astigmatism in the instrumental point spread 
function, the instrumental line width ai ns t r depends 
on the wavelength and on the X spatial direction, 
but not on the Y direction. We constructed maps 
of GinstriXiX) which were subtracted in quadrature 
from the data to obtain the intrinsic velocity disper- 
sion of the gas. At our observed wavelength of Ha 
(~6800A), cjj ns t r .~51~-62kms _1 , being close to the 
maximum in the centre of the field. The total er- 
ror on the velocity dispersion accounts for the uncer- 
tainty of ~5kms _1 in (Ti nstr . The errors from the 
fit on radial velocity are very small (<5kms -1 ), so 



SNR=5 




6750 6800 6850 6900 6950 7000 



wavelength 

Figure 4: Measurement of the emission lines in the 
group [Nn]-Ha-[Nn] for two values of SNR of Ha. 
From top to bottom the real data and one of the 
simulation in a spaxel with SNR~ 5 (first two plots) 
and in a spaxel with SNR~ 50 (last two plots). The 
black curves are the observed and simulated spectra, 
while the red curves are the fits to the emission lines 
and the continuous lines are the fits to the continuum. 
SNR=50 is the mean value of our data. 
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the main contribution comes from the uncertainty in 
wavelength calibration, which is ~13kms _1 . 

As explained above, we chose to measure all lines 
independently in order to be able to assess the reli- 
ability of our flux ratios. With this aim, we com- 
pared our measurements of the [Nil] AA6583/6548 
and [Om] AA4959/ 5007 line ratios to their the oret- 
ical v alues of 2.976 (IDopita fc Sutherland! . l2003h and 
2.936 (jStorey &: Zeippenl . l2000h . respectively. We found 
that for an error on the flux ratio less than 30%, the 
involved lines must have at least SNR>10 (consistent 
with the above errors on the individual fluxes), and 
adopted a more conservative limit of SNR>20. 

To achieve this SNR, we applied a spatial adap- 
tive binning to our data by means of the 'Weighted 
Voron oi Tessellation' (WVT) by Diehl and Statler @ 
( 200 6]) , which is a ge n eraliz ation of the algorithm that 
Cappellari fc Copinl (|2003h developed for SAURON 
data. The WVT performs the partitioning of a re- 
gion based on a set of points called 'generators', the 
points around which the partition of the plane takes 
place. The partitioning is repeated iteratively until 
some target SNR is achieved in all bins. The WVT 
algorithm allows to 'manually' set a number of gen- 
erators, which we chose on the basis of the position 
with respect to the stellar disc. We obtained a total 
of 104 regions. As input data we used the signal and 
noise of the Ha emission line, which has the widest 
spatial coverage. The faintest line involved in our 
flux ratios is the [Oi] A6300 line, whose flux is typi- 
cally about ~ 20% of Ha, so that we set SNR=100 
for Ha as the target SNR for the WVT algorithm in 
order to achieve the required SNR for the [Oi] line. 
The WVT is adopted to derive flux ratios, dust at- 
tenuation, and SFR, but not for the kinematics, for 
which no spatial binning was necessary. 




Figure 
tours 



5: 

of Ha 



5-10 16 ergcm 2 s 1 A 1 arcsec 2 



Logarithmically-spaced con- 
flux in the range 5-10 -18 - 
from WiFeS are 
superimposed on the ESO-WFI 72-band image of 
the galaxy. Extraplanar gas is visible out of the 
disc in the NW (upper right). The scale is shown 
on the bottom. The Ha flux was corrected for dust 
attenuation (Sect. HT3l) . 



6 Results 

Figure [5] shows the contours of the Ha emission-line 
flux from WiFeS tracing the ionized gas superim- 
posed on ESO-WFI i?-band image. The ionized gas 
extends for about 13kpc in projection out of the 
galaxy disc in the NW direction. No gas is detected 

4 http://www.phy.ohiou.edu/diehl/WVT 
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beyond that distance, although pointing #2 covers has an intrinsic diameter of ~0.35" (350 pc). Alter- 
further 13 kpc North from the most external gas isophotffiiatively, SOS 115228 could be a much more distant 



in Fig. [5j In order to understand the origin of the 
extraplanar gas, we examine the structure and mor- 
phology of the ionized gas, derived the gas and stellar 
kinematics and the dust extinction across the galaxy. 

6.1 Morphology of the Ha emission 

The Ha narrow-band image allows us to resolve the 
structure of the extraplanar gas. Fig. [6] shows the 
MMTF Ha image (white) combined with the UKIRT 
if-band (red) image (left panel) and with the con- 
tours of Ha flux measured with WiFeS (right panel) . 
The most noticeable feature of the Ha image are the 
compact (<^0.5"; ^500 pc) knots of Ha emission seen 
all along the NW side of the disc. Some of these 
knots seem "tethered" to the disc by faint filamentary 
strands. The most distant knot is seen at ~ 13 kpc 
(in projection) from the galaxy major axis, but most 
knots are much closer, within 3-4 kpc of the disc. All 
of these knots and filamentary structures are com- 
pletely absent in the MMTF continuum image, con- 
firming that this is from Ha emission with little if 
any underlying stellar continuum component. It is 
notable that there are no such H-alpha features on 
the SE side. 

The most luminous knot has an Ha luminosity 
of 7.2 xlO 38 ergs -1 , while the faintest detected knots 
have L(.ffa)~1.5x 10 38 ergs -1 . If all of this emission 
were powered by star formation, these luminosities 
would correspond to SFRs of 0.0008- 004 M ^ yr" 1 
based on the calibration of Kennicutt ( 19981 ) with a 
Kroupa IMF (see Sect. 17. 1|) . A bright, clumpy Ha- 
emitting region ~6kpc in extent is located in the 
disc ~12kpc SW of the galaxy centre, indicative of a 
localized ongoing starburst (see Sect. I6.3| ). 

The small galaxy SOS 115228 (Fig [H right) SE 
of SOS 114372 has also an Ha-emitting region as- 
sociated with it, which appears ~2.8x brighter in 
our MMTF Ha filter than the continuum filter. Al- 
though SOS 115228 does not have a known redshift, 
this might be suggestive of a starbursting dwarf galaxy 
at the same redshift (within ~400kms _1 ) as 
SOS 114372. Although clearly compact, SOS 115228 
does not appear as a point source in Ha, but likely 



active (starburst or AGN) galaxy whose H/3 emission 
happens to be detected in our MMTF image. The 
redshift of this galaxy would then be z~0.42. This 
difference in redshift would be consistent with the 
~5 mag difference in if-band magnitude with SOS 114372, 
if the two galaxies have comparable masses. Unfor- 
tunately, no conclusion can be drawn on the nature 
of this object without spectroscopic data. 

6.2 Gas and stellar kinematics 

The gas velocity field derived from IFS observations 
is shown in Fig. [7] (left panel). The white contours 
trace the -R-band stellar continuum derived from the 
spectrum once the emission lines have been removed. 
The dashed lines mark the major and minor axes of 
the disc extending to 3xr^, and cross at the if-band 
photometric centre. The black curves are iso-velocity 
contours. The kinematic centre is assumed to co- 
incide with the photometric centre of the if -band 
image. It turns out that this choice also produces 
the smallest asymmetries in the major axis kinematic 
profile. 

The velocity field is complex. Within the disc, the 
overall appearance of the field is that of a rotating 
disc, although significant departures from a simple 
rotation are present everywhere, as it is also clearly 
shown in the radial velocity profiles along the major 
and minor axes (Fig. [8]) . 

The major axis velocity profile (Fig. [5]) is fairly 
symmetric in the inner ~5 arsec. Beyond that ra- 
dius, the radial velocity in the NE side of the disc 
increases until it reaches ~250kms~ 1 at ~16 arcsec 
from the centre and then it starts decreasing until the 
last observed point at 21 arcsec. In the SW side, the 
absolute value of the velocity first decreases to a local 
minimum of ~140kms _1 at 10 arcsec from the cen- 
tre and then slightly increases until the last observed 
point. The whole situation depicted in Fig. [7] is even 
more complex than this, as is immediately clear from 
the shape of the iso-velocity contours. The velocity 
field in the disc is asymmetric also with respect to 
the minor axis, with the SE side of the disc generally 
more red-shifted than the NW side. There are two lo- 
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Figure 6: Left: Composite image of SOS 114372 consisting of the MMTF Ha image (white) over the UKIRT 
iT-band image (red). SOS 115228 is seen as the white spot closest to the bottom. The Ha flux is not 
corrected for dust attenuation. Right: Contours of Ha flux measured with WiFeS over the MMTF Ha 
image. The scale for the two panels is indicated. 



cal maxima in the SE side, one at ~3 NE from the 
centre and one at ~2r^ SW from it. The kinemat- 
ics of the extraplanar gas appears to be remarkably 
dominated by the rotation characterising the disc, 
maintaining the iso-velocity contours continuous up 
to the most external limits. We notice that, NE from 
the galaxy centre, the extraplanar gas is in the av- 
erage blue-shifted with respect to the nearest gas in 
the disc, while in the SW it is red-shifted, thus pro- 
ducing 'fan-shaped' iso-velocity contours. The minor 
axis radial velocity profile allows us to identify a gas 
stream approaching the observer, reaching a velocity 
of ~100 km s -1 at ~3-4 arcsec SE from the nucleus. 

From the major axis radial velocity profile we de- 
rive a dynamical mass within 20 kpc of A / ldj /n =2.1-10 11 

The gas velocity dispersion field is also quite com- 
plex (Fig. right panel) . We first notice that the ve- 
locity dispersion a rarely reaches values as small as 
would be expected in normal turbulent HII regions 
(<7~20-30kms -1 ). The velocity dispersion is partic- 
ularly high (> 100 kms -1 ) over a significant fraction 
of the extraplanar gas. High values of a are also ob- 
served in two regions near the SE border of the disc. 



The most remarkable of the two is the triangular area 
with a ^120 km s -1 extending SE from the nucleus, 
where we also measure the absolute maximum of the 
velocity dispersion (~180kms _1 ). This area corre- 
sponds to the stream directed towards the observer 
noted above. We remark that the value of the veloc- 
ity dispersion also depends on the spatial sampling of 
the spectrograph, increasing with pixel size and see- 
ing, as the combination of the motions of more gas 
elements is measured in each resolution element. This 
is certainly a major origin of our average high values 
of a, because it is clear from the radial velocity field 
that complex motions are taking place giving rise to 
different possible kinds of superimpositions. We no- 
Mjce that our resolution element of 1" x 1" corresponds 
to ~lx7kpc 2 once projected on the plane of the disc. 
We also remark that the presence of large quantities 
of obscuring dust may complicate the interpretation 
of the gas kinematics, since dust may selectively hide 
parts of the mo ving g as (e.g. Giov anelli & Haynes 
20021 : Baes etal. 120031 : Kregel et al. 120041 : Valotto & 



Giovanelli [2004J), although this would probably not 
increase the measured velocity dispersion. 
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Figure 7: Gas kinematics of SOS 114372 derived from the fit to the Ha emission line. The light gray rectangle 
marks the region not covered by the two pointings. The white iso-density contours trace the galaxy red 
continuum and the dashed lines, extending to 3xr^, show the positions of the minor and major axis and 
cross at the -ftT-band photometric centre. Left: Gas velocity field. Each pixel is colour coded according 
to the measured radial velocity relative to the galaxy centre. The black curves are iso-velocity contours. 
Right: Gas velocity dispersion. Each pixel is colour coded according to the measured velocity dispersion. In 
this figure and in the following maps the ticks on the axes indicate the apparent distance from the -fT-band 
photometric centre in arcsecs. The velocity scales are on the right of each panel. 
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Figure 8: Gas line-of-sight velocity (upper panels) and velocity dispersion (lower panels) profiles along the 
major (left panels) and minor axis (right panels). The red lines represent the same data reflected with 
respect to the origin, to facilitate the inspection of asymmetries. Distances from the photometric centre are 
positive towards NE (major axis) and NW for the minor axis. Directions in the sky are given on the top of 
the figure. 
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Figure 9: Top: The stellar velocity field with Iso- 
velocity contours. Bottom: The stellar radial velocity 
profile along the major axis (red dots) with its reflec- 
tion with respect to the centre superimposed (green 
curve) . 



10 




-15 -10 -5 



0.0 0.4 0.8 1.2 1.6 2.0 2.4 2.8 3.2 



Figure 10: Dust attenuation derived from the Ha/H/3 
line ratio. 

From the fit of the stellar continuum we derived 
the velocity field of the stars. In order to estimate the 
uncertainties in the stellar velocity we used simula- 
tions in which random noise is added to the spectra. 
A typical uncertainty of ~2kms -1 for spaxels in the 
disc is found. This low uncertainty is probably due to 
the previous 3x3 smoothing in the spatial direction. 
Including the uncertainty in wavelength calibration, 
this sums up to ~14kms _1 . 

The derived stellar radial velocity field and major 
axis radial velocity profile are shown in Fig. [9l The 
iso- velocity contours of the stellar component appear 
symmetric with respect to the centre, while they are 
skewed with respect to the minor and major axis. 
Along the major axis, the radial velocity grows lin- 
early and rapidly in the inner ~1. 5-2 arcsec, where 
there is an abrupt change of slope followed by an- 
other rapid increase. At ~5 arcsec from the centre the 
radial velocity assumes an almost constant value of 
~180kms -1 until the last observed points. The kine- 
matic centre of the stellar velocity field is displaced 
0.3 arcsec (280 pc) NE with respect to the photomet- 
ric (K band) centre of the galaxy (and consequently 
to the kinematic centre of the gas). However, we will 
not further discuss this point, since this displacement 
is smaller than our resolution element. 
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6.3 Dust extinction across the galaxy 

Figure [TU] shows the distribution of the dust atten- 
uation derived from the Ha/H/3 line ratio, in terms 
of the visual extincti on Ay. We use the theor etical 
attenuation curve bv lFischera Dopital (|2005h with 
i?y=4.5. The dust extinction is estimated in each of 
the binned region of Fig. [TT1 The level of dust attenu- 
ation is notably asymmetric along the major axis, be- 
ing universally lower on the NE side of the disc than 
the SW side. The highest dust extinction (Ay~2.4- 
3.2 mag) is observed in the centre of the galaxy, ex- 
tending along the minor axis to the SE edge. Highly 
attenuated regions (Ay~1.9-2.6 mag) are also found 
in the SW disc ~12kpc from the galaxy centre and 
in the extraplanar gas. This extinction map agrees 
with what can be inferred from the optical images, 
where dust absorption is clear in the galaxy centre 
(see also Fig. [1] right panel). 

7 Properties of the gas and star 
formation across the galaxy 

Our data, complemented by available multi-band imag- 
ing (from far-UV to far-IR), allow us to investigate 
the nature of the gas and the star formation across 
the galaxy, but first we examine the gas and stellar 
kinematics. 

The stellar kinematics is typical of a barred galaxy, 
with the twisting of the iso- velocity contou rs due to 



also for much more severe binning (up to 7x7 pixels). 
The complexity of the gas velocity field cannot be ex- 
plained by the presence of the bar alone, as can be 
seen by comparing the iso- velocity contours in Fig. [7] 
to the veloci ty fields of bar r ed galaxies repor t ed, fo r 



instance, in 



Bosmal dl98lh. ICourteau et al.1 faooi ) 



Fathi et al.1 (j2005h and Ide Narav et HI (fcOQflh . All 



these considerations make us to conclude that the 
gas kinematics is heavily perturbed, while the stellar 
kinematics is not. 

7.1 Physical properties of the gas 

The key diagnostics for the examination of the mode 
of excitation of an ionized pla sma were introduced by 
Veilleux fc Osterbrockl (jl987h . These use the [Nil]/ 
Ha, [Sii]/Haor [Ol] /Ha ratios plotted against the 
[Oiii]/H/3 ratio, and have the great advantage of be- 
ing affected very little by dust extinction. They allow 
a detailed classification of the excitation as either due 
to young stars, or due to an active nucleus. The clas- 
sification scheme has been refine d by Kewley et al. 
), Kauffmann et al. (|2003l ) Kewley et al. 



2001 



the n on-circular motions induced by the bar (lAthanasso 
20021) and with a 'double-bump' rota tion curve ( Chung 



2004; iBureau h Athanassoulal . 20051 ) . A detailed anal- 
ysis of the stellar kinematics is beyond our goals, 
since it would requir e muc h higher spatial resolution 
(e.g. Emsellem et al. 20071 ). but the important point 
for this work is that the stellar velocity field appears 
to be regular and symmetric, without any signs of 
perturbation. One may ask whether the regularity 
of the stellar kinematics field, in sharp contrast with 
the gas, is due to the smoothing applied to the spec- 
tra to increase the SNR of the stellar component. We 
verified that this is not the case by applying the same 
3x3 binning to the gas velocity field. The effect was 
just a smoothing of the small-scale features, with the 
gross properties remaining unchanged. This holds 



(|2006l ) by the use of more sophisticated models to 
define the allowed range of ratios produced in HII re- 
gions, and using the SDSS galaxy sample to create 
a semi-empirical classification to distinguish between 
Seyferts or LINEFi-like AGN. 

We investigate the physical properties of the gas 
using the pointing #1 where the suitable SNR (SNR=100 
hjr ffifl/fil^p ^ • ^ s achieved by binning the data 
Liiujugn une vv \fT. In Fig. [TT]we show the 104 differ- 
enFgafaxy regions where independent measurements 
of the spectral indices were derived. These include 
regions outside of the disc where the extraplanar gas 
is detected although with a lower SNR per spaxel. 

In Fig.[12]we show the measured line ratios for the 
different galaxy regions, colour coded as in Fig. [TTJ 
along with the theoretical predictions and limits for 
photoionization and shock excited models. The red 
curves mark the theoretical upper limit for HII re- 
gions, while the dashed blue curve in the left panel 
delineates the empirical limi t to st ar-forming regions 
defined by Kauffmann et al. ( 20031 ). Between the red 
curve and the dashed blue curve, 'composite' systems 
are found whose emission is characterized by a mix- 



20 



1.5 i i i i i | i i i i | i i i i | i i i i | i i i i | i i i i i i | i i i i | i i i i | i i i i | i i i i | i i i i | i i i i | i i i i | i i i i 



1.0 - 




log([NII]/Ha) log([SII]/H a) log([OI]/H a) 



Figure 12: Line flux diagno stic diagrams of t he different regions of SOS 114372 (see text). The shock and 
photoionization models by lEich et alJ (j201lh are superimposed on the measured flux ratios. The black 
curves are drawn for different fractions (from to 1) of Ha flux contributed by shocks as indicated. For 
comparison the theoretical (red curve) and empirical (blue dashed curve in the left panel) upper limits for 
HII regions are also indicated as well as the separation between AGN and LINER (solid blue line in the 
central and right panel). 




Figure 11: The 104 galaxy regions identified by the 
WVT are shown in different colours. For each of these 
regions reliable and robust measurements of emission 
line ratios are derived. The colours used in this figure 
are consistently used in Figs. [T2land[T71 



ture of AGN and HII regions. The blue continuous 
line in the central and right panels divides the line 
ratios typical of Seyfert galaxies (upper part of the 
diagram), from those of low ionization emission line 
regions ('LINER'). These curves are extensively used 
in the literature and plotted here to facilitate com- 
parison with other works. 

The black curves represent shock a nd photoion- 
ization models from Rich et al. ( 201ll ). in which a 
fraction of the radiation (as indicated by the num- 
bers close to the curves) is produced by shock excita- 
tion and the remaining comes from photo-ionization. 
These models are computed using the MAPPINGS 
Illr code, t he latest version of the code originally in- 
troduced in lSutherland Dopita (Il993h. The abun- 



dance set is taken from Grevesse et al 



jpOld ) and 



wc 



use t he standard (solar reg ion) dust depletion factors 



from iKimura et al.l (|2003l ) . The shock models have 
velocities in the range 100 to 200kms _1 in steps of 
20kms _1 , with hydrogen assumed to be fully pre- 
ionized. In using these depletion factors we are im- 
plicitly assuming that the shock is not sufficiently fast 
to sputter the dust grains advected into the shocked 
region. The models presented here have Z/Zq=2 
(log(0/H)+12= 8.99), and have a transverse mag- 
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netic field consistent with equipartition of thermal 
and magnetic energy (B=5 /j,G for n c =10 cm~ 3 ; Bocn c 1 

The photoionization models are generated using 
UV flux distributions appropriate to continuous star 
form ation computed using the Starburst99 code (jLeithe 
1999 1) . We have use d the older models employed by 
Kewlev et all (|200lh rather than more recent ones, 
since the hardness of the radiation field provides a 
closer match to the HII region line ratios observed in 
the SDSS galaxy sample. We have also adjusted the 
nitrogen abundance to provide a still closer match to 
the SDSS galaxies. The photoionization models are 
characterized by an ionization parameter 6.5< log[g]<7. 
This ionization parameter can be converted to the 
more commonly used dimensionless ionization param- 
eter U by U=q/c, where c is the speed of light. 

We have created mixing models between each pho- 
toionization model, and all shock models. These are 
characterized by a fraction / of the flux at Ha being 
contributed by shocks. These models do not create a 
simple mixing line on the diagnostic plots, but rather 
a zone within which the line ratios are consistent 
with a certain mixing fraction of shock excitation. 
In Fig. [12] we show the average of the line ratios pro- 
duced by mixing fractions, /=0, 0.15, 0.4, 0.7 and 
1.0. Note that the line ratios are sensitive to quite 
small contributions from shock emission - as little as 
15% contribution by shocks can change observed line 
ratios by a factor of two. 

The observed line ratios cannot be simply ex- 
plained by HH-like photoionization, but display a 
spread which is characteristic of a mixture of both 
shock excitation and photoionization. All parts of 
the galaxy seem to be contaminated by at least a 
small fraction of shock excitation, 0.05^/^0.1. The 
regions which are dominated by star formation and 
the associated HII regions are preferentially located 
in the SW portion of the galaxy (purple points). Re- 
gions in the extraplanar gas to the NW of the galaxy 
(red/orange points) have much higher shock excita- 
tion fractions, ranging up to 0.4-0.8. In addition, 
there is a strong concentration of shock-excited gas 
along the minor axis SE from the nucleus (light blue 
points) and in another Northern region of the SE edge 
(green point). 
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Figure 13: The logarithmic SFR density across the 
galaxy derived from Ha flux corrected for the ex- 
tinction and shock-excitation, where necessary. The 
contours show the isophotes of the 24/im flux. The 
three boxes outlined by white dashed sides mark the 
regions of the galaxy where the stellar population 
analysis was performed. The three regions are, from 
right-bottom to top- left: the starbursting region in 
the SW disc, the central bulge region, and the NE 
disc. 

7.2 Star formation across the galaxy 

To study the star formation across SOS 114372, we 
first analyze the ongoing star formation from the flux 
of the Ha emission line. Then we study the recent 
star formation history from the distribution of the 
ages of stellar populations in three distinct regions of 
the galaxy - a starbursting region in the SW disc (see 
below), the central bulge region, and the NE disc (see 
FigUBJ. 

7.2.1 Ongoing star formation 

We derive the current SFR from the Ha flux tak- 
ing int o account t he eff ects of dust extinction fol- 
lowing Kennicutt (1998). For the derivation of the 
SFR we approximately removed the contribution of 
the shock-ionized gas to the Ha flux - the Ha fluxes 
of the galaxy regions with >50% of emission coming 
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from shock-ionized gas are reduced by a factor of two. 
This correction is however very small compared to the 
uncertainties. The logarithm of SFR surface density 
across the galaxy is given in Fig. [TBI Intense star 
formation is found in the disc spanning from 0.01 to 
0.34 Mq yr _1 arcsec -2 with the global maximum lo- 
cated in the galaxy centre. There is also a notable 
local maximum in the SW region of the disc at ~ 12 
kpc from the centre, with SFR=0.2 M Q yr~ 1 arcsec -2 . 
The contours in Fig. [13] show the distribution of the 
24/xm flux which, in agreement with the SFR derived 
form Ha, presents the maximum in the centre and the 
second peak in the SW side. These two regions are 
also the most obscured regions identified in Fig. [10J 
These results imply a coherent scenario where the 
current SFR, derived from the Ha emission and the 
24/iin flux, occurs in the highly dust-extinguished 
molecular clouds. The integrated Ha-derived SFR 
of SOS 114372 (adopting the Kroupa IMF) amounts 
to 7.2±2.2 Mq yr _1 . The error takes into account 
the uncertainties related to the flux and attenua- 
tion measurements added to a 30% uncertainty due 
to different calib rations of Ha as a SFR indicator 
( Kennicutt . 19981 ). From the difference between raw 
and attenuation-corrected Ha fluxes we can estimate 
the obscured SFR as 5.3±1.6 M yr _1 . This obscured 
SFR is also measured by the IR dust emission turning 
out to be 6.65^^19 M yr _1 . The two independent es- 
timates of the obscured SFR are therefore fully con- 
sistent. 

7.2.2 Recent star formation history 

Having seen that the distribution of dust obscuration 
and ongoing star formation is strongly asymmetric, 
being significantly reduced in the NE side of the disc, 
while the SW side shows an ongoing highly obscured 
starburst, we now examine the three regions of the 
galaxy indicated in Fig. [13] (SW disk starburst, cen- 
tral bulge starburst, NE disk) separately. To deter- 
mine the ages of stellar populations, we fit the spec- 
trum with a non-negative linear combination of 40 
SSPs, this time with attenuation of the stellar com- 
ponent by dust also left as a single free parameter (we 
do not consider selective attenuation for the youngest 
stellar populations), resulting in the best-fit model 



for the stellar spectrum. This is then subtracted, 
leaving the residual noise and emission spectrum. To 
achieve the signal-to-noise ratio needed for this anal- 
ysis, we added the spectra within each of the three 
regions indicated in FigJTBI 

In the SW disc starburst region (Fig. [T4"|) . the 
shape of the continuum, together with the weak 4000A 
break (d n 4000=1.09), indicate that the spectrum is 
dominated by young stars. This is confirmed if we 
split our best-fit model into its contribution from 
young (<lGyr; thin blue curve) and old (>lGyr; 
thin red curve) stars. Young stars contribute ~60% 
of the blue luminosity and ~3.5% of the stellar mass 
in this region. Of this young (<lGyr) component, 
35% of the stellar mass and 60% of the luminosity 
comes from stars younger than ~100Myr. This im- 
plies a starburst characterised by a ~ 5x increase 
in the SFR over the last ~100Myr. This burst is 
ongoing as implied by the very strong (despite be- 
ing highly obscured with Ay~1.9-2.6; Fig. 10) Ha 
emission in this region with EW(Ha)~120A. Further 
Balmer emission lines are visible in Fig. [Til from H/3 
all the way up to H£. The old stellar component (red 
curve) can be simply characterized as a single 10-15 
Gyr old population. 

In the central bulge region (see Fig. (2{H, the shape 
of the continuum appears dominated by an old stellar 
population. The splitting of our model stellar pop- 
ulations into old and young, shows that young stars 
contribute just ~30% of the luminosity and 1.0% of 
the stellar mass. Moreover, young stars are predomi- 
nately of age ^100 Myr, as also suggested by the over- 
all smooth appearance of the continuum from this 
component. 

In the NE disc region (Fig. fT5j) . the stellar con- 
tinuum is dominated by young (<lGyr) stars, con- 
tributing to ~60% of the luminosity and ~4% of the 
total stellar mass. Our stellar population modeling of 
the young component suggests multiple populations 
spread over the full range 60-1000 Myr with just ~5% 
of the stellar mass in this young component coming 
from stars younger than 100 Myr. This suggests that 
the average SFR over the last ~100 Myr is ~2x lower 
than that averaged over the last Gyr. 



5 Since the results of the present fit are qualitatively similar 
to those shown in Fig. [5] we refer to that figure in this analysis. 
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Figure 14: The integrated spectrum (black curve) and full population synthesis fit (magenta) for the starburst 
region in the SW side of the disc. The contributions to the model fit from the old (>1 Gyr; thin red curve) 
and young (<lGyr; thin blue) stellar populations are shown. The shaded regions indicate the wavelength 
ranges excluded from the fitting process. The residual emission component after subtraction of the model 
stellar continuum is shown by the thick blue curve. 
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Figure 15: The integrated spectrum (black curve) and full population synthesis fit (magenta) for the NE 
side of the disc. The curves are coloured as in Fig 1141 The line indices discussed in the text are indicated. 
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To better understand this situation, we compare 
the stregths of [On] and H/3 emission lines, which 
trace the ongoing star formation, to the H<5 absorp- 
tion line, tracing the intermediate-age (~1 Gyr) A- 
type stars. We measure EW([Oll])=-5.2±0.06A, 
EW(H(5)=5.24±0.08A and EW(H/3)=-2.19±0.06A. 
This combination of line strenghts places the NE disc 
on the boundary between the weak post-starburst 



accounts for ~20% of the total SFR, our data im- 
ply a starburst characterised by a 5x increase in the 
SFR over the last ~100Myr. This burst is ongo- 
ing. In the NE disc, our data still show ongoing 
star formation, but significantly lower than in the 
rest of the disc. The full spectral modeling and the 
line strengths show that we are observing a 0.3 Gyr 
starburst immediately after it has been shut down. 



(k+a) and e(a) (post-starbust with residual star-formation) 
classes 



(see Dressier et al. . 19991 ). The mean relation be- 



8 Origin of the extraplanar gas 



tween EW(H5) and EW([Oll]) derived by |Dressler et al. T he mogt prominent and notable physica i character 
(|2004l ) for normal star-forming galaxies predicts an 
EW(H(5)~2.0A| for a galaxy with EW([0n])=-5 A. 
The fact that our observed level of H<5 absorption is 



somewhat higher than expected again indicates that 
the current SFR is significantly lower than that av- 
eraged over the previous Gyr. 

As a further probe of the ong oing-to-recen t SFRs 



we derived the Rose Can index (jRosd . 1 19851 ) , which 
is defined as the ratio of the counts in the b o ttoms 
of the Can H+He and Can K lines. iRosel (j 1985h 
showed that the Call index is constant at ~1.2 for 
stars later than about F2, but decreases to a mini- 
mum of ~0.4 for AO stars as the Can lines weaken, 



befor e increasing again as the He line fades. Further- j gg Cortese et al 



Leonardi fc Rose! (|l996l l combined the Can in 



more, 

dex to the H5/Fe I A4045 index to determine both the 
duration and the strength of starbursts. We measure 
a Can index of 0.635 and a H5/FeiA4045 index of 
0.745. These values confirm that our spectrum is 
dominated by A-type stars and, more importantly, 
show that we are observing a starburst which lasted 
for ~0.3 Gyr and which has been just shut down (cf. 
Fig. 4 of iLeonardi fc Rose! fl99fih . 

In summary, the total SFR from the Ha flux 
amounts to 7.2±2.2 M Q yr _1 fully consistent with the 
previous estimates from UV and FIR measurements. 
The star formation across the galaxy mainly takes 
place in two regions - in the central bulge region and 
~ 12kpc SW from the centre. In the central bulge 
region there is an intense ongoing, heavily obscured 
star formation accounting for ~30% of the total SFR 
of SOS 114372. In the SW starbursting region, which 



istic of SOS 114372 is the extraplanar ionized gas ex- 
tending out to a projected distance of 13kpc NW 
from the disc over its full extent following the ro- 
tation of the disc. This gas is resolved into a com- 
plex of compact Ha-emitting knots and faint filamen- 
tary structures and shows a high fraction (0.4-0.8) 
of shock excitation. All this supports a scenario in 
which gas is being driven out from the galaxy disc. 
The high levels of dust attenuation with Ay ~ 1.9 mag 
also found associated with the extraplanar gas, sug- 
gests that the dust is stripped and dragged away from 
the disc together with gas, in agreement with the re- 
duced _djj^t_djs^_obs^rved_in Hi-deficient Virgo galax- 
" (|2O10h . 



6 with 95% of normal star-forming galaxies lying in the range 
0.0-4.0 A 



Causes of large-scale outflows of gas from a galac- 
tic disc may be galactic winds produced by AGN 
or starbursts, tidal interaction or ram-pressure strip- 
ping. In the following we will outline the arguments 
that lead to the idea that RPS is actually the origin 
of the observed gas outflow. 

8.1 Galactic winds 

We notice that the outflow takes place almost uni- 
formly all over the disc and is observed only in one 
direction (NW, see Figs. [5] and [6|). Since it extends 
in projection several kiloparsecs outside of the disc, 
dust attenuation could not hide a possible (symmet- 
ric) component on the SE side. The observed outflow 
is therefore intrinsically asymmetric. Galactic winds, 
instead, produce bipolar outflows, with a structure 
ranging from biconical to egg-shaped, originating from 
a nuclear starbu rst or AGN, and ext ending on both 
sides of the disc (|Veilleux et al.l . 120051 . and references 
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therein). We therefore exclude a galactic wind as the 
cause of the bulk of the gas outflow observed here. 

Nevertheless, we don't exclude that a galactic wind 
could actually be in place close to the nucleus, where 
star formation is taking place at a rate of ~2 M Q yr _1 
and the Ha image shows hints of gas escaping from 
the centre of the SE disc (see Fig. EJ). Let us con- 
sider the with very high values of the velocity 
dispersion 120 km s _1 , Fig. [7]) located along the 
minor axis SE from the galaxy centre. We notice 
that this area has a roughly triangular shape, with 
a vertex near the galaxy centre. The minor axis ra- 
dial velocity profile (Fig. [8]) shows that the gas is 
moving towards the observer, with a speed steadily 
rising up to ^lOOkms -1 at ~ 3kpc from the centre. 
The high values of a suggest the presence of different 
components of gas with different radial velocities, so 
that the measured ^lOOkms" 1 is an average velocity 
(weighted by the H Q flux). 

All these features, together with the fact that the 
radiation here is shock-dominated (Sect. 17. ip . sug- 
gests that the triangular region SE from the centre 
may be a semi-cone prod uced by the nuclear star- 
burst (jVeilleux et all 120051 ) . The absence of the com- 
panion semi-cone in the NW direction could be ex- 
plained by dust absorption if the disc is oriented with 
the SE side furthest from us. We will add further 
support for this orientation in Sect. [9j Given the in- 
clination of the galaxy and the apparent extension 
of the wind cone of 6.5arcsec from the centre, and 
allowing for large uncertainties due to poor spatial 
resolution, we infer that the cone extends ~ 6kpc 
above the disc, within the typical sizes of observed 



galactic winds (jVeilleux et all 120051 ) . 



As already noted, the area SE from the nucleus is 
the most dust-obscured (excluding the nuclear region 
itself; Fig. flOl) . This would be easily explained by a 
large amount of dust collected in the starbursting re- 
gion and entrained in the wind, which i s a rather com 



mon circumstance in galactic winds (jHoopes et al. 
20051 . and references therein). 



8.2 Tidal interaction 

The best candidates for a major tidal interaction 
are the cluster members SOS 114493 and SOS 112392 
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Figure 16: Close-up of SOS 114493 (left) and 
SOS 112392 (right) from a ESO-2.2m WFI deep ri- 
band image. The intensities are stretched by a 
logarithmic scaling to emphasize the lowest surface 
brightness regions. The deep-red ellipses are drawn 
from the models derived from a simultaneous fit to 
the two galaxies using two Sersic laws . The fit was 
performed with GALFIT (jPeng et all l2O10h . 



The 

surface brightness of both galaxies smoothly follow 
the models without any appreciable sign of perturba- 
tion. One side of this figure corresponds to ~45kpc. 

(Fig. Q] and [To]) with redshifts z = 0.050 and z = 
0.046 and i^T-band magnitudes 13.49 mag and 12.61 mag, 
respectively. These two galaxies are at projected dis- 
tance from SOS 114372 of ~45kpc and ~55kpc re- 
spectively, and their stellar masses are .Mii4493=2-1O 1O M 
and 7Wii2392=7-lO lo M (Paper I). We observe no tidal 
bridge connecting SOS 114372 to either of these ob- 
jects nor any signs of disturbance in their i?-band 
images, as is shown in Fig. [TBI The direction of 
the gas outflow from SOS 114372 does not appear 
to be related to any of these galaxies. Furthermore, 
no obvious sign of perturbation is present in any of 
SSC galaxies within 200 kpc projected radius from 
SOS 114372. 

As already noted in Sect. 16.21 the stellar veloc- 
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ity field of SOS 114372 displays no sign of distur- 
bances, indicating that the galaxy is not suffering 
any significant external gravitational perturbation. 
The lack of a tidal tail, the undisturbed stellar disc, 
the existence itself of the gas disc, and the absence 
of perturbed neighbours make us exclude the occur- 
rence of tidal interaction, unless it is at a very early 
stage, with the closer galaxy (at a distance anyway 
> 45kpc) just starting its first approach. N umeri- 
cal simulations ( Kronberger et all 2006 . 20071 ). how- 
ever, do not predict gas outflow at this stage, either 
for two galaxies of the same mass (as in the case 
of SOS 114372-SOS 112392) or for the unequal-mass 
case (as for SOS 114372-SOS 114493). 

The next candidate for a tidal interaction is the 
faint object located SW of SOS 114372 (SOS 115228, 
see Fig. [1] right panel). Although we do not have 
redshift information for this galaxy, it appears much 
brighter in our narrow-band Ha image than in the 
continuum, which would suggest that it is a star- 
bursting dwarf galaxy at a redshift within a few hun- 
dred kms -1 of SOS 114372. As tidal torques are 
effective in channelling rapid gas inflows capable of 
fuelling nuclear starbursts, it seems reasonable that 
an interaction between SOS 114372 and SOS 115228 
could be behind the observed starburst in the dwarf 
galaxy. It is much less likely that this interaction 
could induce the large-scale gas outflow in SOS 114372 
however. Firstly, the .fT-band magnitudes of the two 
galaxies differ by ~ 5 mag, and we found a factor 
about 360 difference in stellar masses (Paper I). This 
makes it unlikely that this object can have such a 
large-scale effect on the gas of its much more massive 
neighbour. Secondly, simulations of unequal mass 
mergers predict gas streams connecti ng the two galax- 
ies (see Fig. 6 of Kronberger et al. 20061 ). while we 
observe the gas outflow in the opposite direction. We 
find no evidence at all for a tidal bridge connect - 
ing the two galaxies, although iKapferer et al.l (|2009al ) 
find that in clusters the impact of pressure from the 
ambient ICM can destroy these tidal bridges, as well 
as further enhance the levels of star formation in the 
interacting galaxies. Notice, however, that the simu - 
lated galaxies considered by Kronberger et al. ( 20061 ) 
to investigate the unequal mass merger have at most 
a mass ratio equal to 8:1. It may be that the notable 



starburst on the SW-edge of SOS 114372 is caused 
by such a combination of the tidal interaction and 
ram pressure stripping. Finally, as already noted in 
Sect. 16.11 it is possible that SOS 115228 is at a com- 
pletely different redshift from SOS 114372. 

All the above considerations and most importantly 
the symmetric stellar distribution of the three brighter 
galaxies lead us to conclude that tidal interaction is 
not the cause of the observed gas outflow. 

8.3 Ram-pressure stripping 

The large-scale gas outflow detected by the IFS data 
is resolved in the narrow-band Ha imaging into a 
complex of compact knots of Ha emission, often at- 
tached to faint filamentary strands which extend back 
to the galactic disc (Fig. [6]). Similar complexes of Ha 
or UV-emitting knots and fila ments have been seen 



for numerous other galaxies (e.g . iGavazzi et all 12001 



Cortese et all. l2007h lYagi et all Eoid : IVollmer et al 



20ld : lFossati et all 12012 )7 The most notable aspect 



of these systems are that all the affected galaxies are 
located in cluster environments, and that the outflow 
is one-sided, often characterized as a cometary tail 
of ionized gas directed away from the cluster centre. 
These two aspects both point to a cluster-specific pro- 
cess in which gas is stripped from an infalling galaxy 
by its passage through the ICM. 

The Ha velocity field of SOS 114372 can be com- 
pared with those of nearby galaxies affected by RPS 



such as NGC 4438 ( Chemin et al 



2005; 



|2009j) and NGC 4522 (IVollmer et al 



200' 



Vollmer et al. 
3) in the Virgo 



cluster. Both galaxies show one-sided extraplanar gas 
following the rotation of the disc, as for SOS 114372. 
The velocity field of NGC 4438 is characterized by 
filamentary structures extending out of the disc. If 
present, such filaments could not be resolved at the 
higher redshift of SOS 114372 with our IFS data. Re- 
markable analogies are found in the Ha images of the 
two g alaxies (compare Fig. 2 of Kenney &: Koqpmann 
1999 with our Fig. [6]). iKennev Koopmannl (|l999l j 
observed that the Ha emission arising from both HII 
regions and diffuse emission is organized in filaments 
extending more than 3 kpc from the outer edge. The 
difference between the two galaxies is in the trun- 
cation radius - NGC 4522 has a truncated gas disc 



27 



at radius of 3kpc, while in the case of SOS 114372 
the ionized gas extends up to 13kpc. This differ- 
ence can be related to different 'ages' of RPS, dif- 
ferent orbits, ICM densities and ICM wind angles. 
We also notice that there are intrinsic differences be- 
tween the two Virgo galaxies and that studied in this 
article. SOS 114372 has an absolute S-band mag- 
nitude of Mb = —21.48, significantly brighter than 
M B = -20.04 and Mb = -18.10 of NGC4438 and 
NGC 4522, respectively (assuming 16 Mpc as distance 
of the Virgo cluster and the magnitudes from Chung 

we are comparing 



et al. 200!) i. This shows that 



galaxies of different masses. Furthermore, NGC 4438 
is located close to the Virgo cluster centre (280 kpc 
from M87, r ~0.3r ?nn) while NGC 4522 is at r ~ r 20 o 
rtUrban e all . l201lh and SOS 114372 is at r ~lMpc 
(~0.5r2oo) from the centre of the rich cluster A 3558. 
It is therefore likely that these galaxies are expe- 
riencing different ICM-ISM interactions. The com- 
mon point between the velocity fields of NGC 4522, 
NGC 4438, and SOS 114372 is that the extraplanar 
regions are still dominated by rotation. 

Hydrodynamical simulations show that ram-pressur< 
and viscous stripping of spiral galaxies produce nar- 
row wakes of Hi gas th at can extend up to 100 - 
200 kpc and last 500 Myr (|Tonnesen fc Brvanl . lioioh . 
These simulations predict a wide range of densities 
and temperatures within these tails, leading to the 
in-situ formation of pressure-supported clouds dense 
enough to form stars. This star-formation plus lo- 
calized heating via compression and turbulence can 
light up the Hi tail with Ha emission, with in tensities 
and morphologies (jTonnesen &: Brvanl. 120121) similar 
to tho se observed in SOS 114372. iTonnesen Bryan 
(l2012l l note that this star formation in the tail does 
not occur due to molecular clouds that were formed 
within the galaxy itself and stripped wholesale, but 
rather because the relatively low-density gas that has 
been stripped can cool and condense into dense clouds 
in the turbulent wake. They also indicate that higher 
ICM pressures create more dense clouds in the wakes, 
resulting in higher SFRs. 

Nevertheless, star formation cannot be the only 
origin of the Ha emission coming from the outflowing 
gas. In Fig. [T7] (upper panel) we focus again on the 
right panel of Fig|12| which plots the [Oiii]/H/3 versus 
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Figure 17: Upper panel: The flux ratios in the dif- 
ferent galaxy regions as in Fig. [T2l The black curve 
corresponds to the shock and photoionization model 
where 50% of the flux is due to shock-ionized gas. 
Central panel: The galaxy regions dominated by the 
contribution of shock-ionized gas colour-coded as in 
Fig. [TTl Lower panel: Normalized distribution of the 
velocity dispersion values. The white coloured his- 
togram refers to the regions dominated by the contri- 
bution of shock- ionized gas, while the gray histogram 
refers to the "low" shock-ionized gas region. The re- 
gions characterized by shock-like line ratios tend to 
2ghave higher velocity dispersions. 



[Oi]/Hq diagnostic. There are seven galaxy regions 
where shock excitation contributes more than 50% 
of the emission. Their flux ratios are located to the 
right side of the black curve in upper panel of Fig. [17] 
and are coloured in the central panel. These regions 
correspond to the gas outflow and to three regions in 
the SE edge of the galaxy. In all these regions, ex- 
cept the easternmost one (shown in light green), we 
also measure strong dust extinction. The normalized 
distribution of the velocity dispersion measured for 
these regions (white coloured histogram in Fig. fT7|) 
is shifted to higher values with respect to the overall 
distribution. We expect the velocity dispersion pro- 
duced by the shocks to be of the same order as the 
shock velocity. If ujcm and tijsm are the number 
densities of the ICM and ISM respectively and if v 
is the relative velocity of the two gas components, 
then the shock velocity is v s =-^JnwM /uism) x v o- 
If nicNi/nishi ~ 10~ 2 -10~ 3 and v Q ~ 1000 km s _1 
then v s ~30-100km s^ 1 , which is what we observe. 
The high velocity dispersions of some non-shocked 
regions may be well explained by the complicated ra- 
dial motions produced by ram pressure and by the 
fact that our resolution element intercepts gas clouds 
for depths of ^ 1/cos i ~ 7 kpc. 

From the distribution of the regions with higher 
velocity dispersion and shock-ionized gas across the 
galaxy, we can infer that SOS 114372 is affected by 
ram pressure which compresses and strips the gas 
out of the galaxy, forming a tail of turbulent shock- 
ionized gas and dust. Part of this gas can cool and 
condense into dense clouds observed in Ha. 

The influence of RPS on the internal kinematics 
has been inves tigated through N-body/h ydrodynamical 
simulations by Kronberger et al. (j2008bl ) demonstrat- 
ing that the gas velocity field shows clear evidence of 
disturbances after a few tens of Myr from the onset of 
RPS. In particular, the rotation curve can be charac- 
terized by asymmetry in the external gas disc while 
the inner disc remains undisturbed. A measure of the 
disturbance is given by the asym metry parameter o f 
the rotation curve as defined by bale et alJ fcOQll ). 
We estimate an asymmetry parameter of about 20% 
for the SOS 114372 major axis radi al velocity profile 



profile of SOS 144372 (extending to ~3.8xr d ; Fig. ED 
is very similar to the gas rotation curve derived for 
the model galaxies after 50 Myr of ram p ressure act- 
ing face-on (Fig. 2 of Kronberger et al. 2008bl ) ex- 
tending to 4xr^. 

SOS 114372 is amongst the brightest SSC galax- 
ies in the mid-infrared with one of the highest SFRs 
measured for the whole sample (see Paper II). The 
galaxy centre and the SW region of the disc con- 
tribute ~30% and ~20% respectively to the total 
SFR measured by the Ha luminosity (i.e. current 
SFR). The location of the star-forming regions in 
the galaxy agrees with the predictions of the sim- 
ulations of RPS where the newly formed stars are 
preferentially located in the compressed central re- 
gion of the galaxy and /or in t h e strip ped gas behind 



the galaxy (IKronberger et all l2008al) depending o n 



the surrounding gas de nsity ( Kapferer et al. . 2009bl ) 



consist ent with what estimated by IKronberger et al 



(|2008bh for RPS. Moreover the major axis velocity 



Sanderson & Ponman ( 2010l ) estimated an ICM den- 
sity of picm = 1-3 x 10~ 28 gcm -3 at the projected 
distance of the galaxy from the cluster core. This 
density implies that all the newly formed stars are 
confined into the disc either because of the low ICM 
density or b ecause the RPS has s tarted less than 
100 Myr ago dKapferer et all l2009bh . The starburst 
detected in the SW disc region produces a ~5x in- 
crease in the SFR over the last ~10 8 yr and could 
be ram-pressure induced. Also in the central bulge 
region, whose spectrum is clearly dominated by old 
(>1 Gyr) stars, we measure a non-negligible contribu- 
tion from young (<^10 8 yr) stars, but this star-forming 
region is most probably related to the presence of the 
bar (see Sect. flOl) . 

The asymmetry observed in the SFR and stellar 
population properties across the galaxy disc can be 
explained in a RPS scenario where the ram pressure 
contributes to the star formation quenching via the 
impoverishment of the gas reservoir, while inducing 
the ongoing starburst in the highly obscured SW side. 
Different stripping efficiencies depending on the sur- 
face density of the stripp ed gas clumps are also found 
by Vollmer et al. ifcood ). 

We can thus conclude that all the observed prop- 
erties of SOS 114372 can be accommodated in the 
RPS scenario. 
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8.4 Environmental conditions for 
ram-pressure stripping 

Now we consider the ICM properties in order to check 

if RPS can be effective for this galaxy. 

For A 3558, Sanderson & Ponman ( 2010l ) derived 
the dark matter density, gas density {picm) and the 
gas te mperature profiles by fitting Ascasibar & Diego 
( 20081 ) models to the Chandra X-ray data. In these 
model s, the dark matter is described by a Hernquist 
(Il990h density profile, which provides simple analytic 
expressions for the mass profile, gravitational poten- 
tial and escape velocity. From this, we can derive the 
ram pressure, P r am=PiCM'v 2 or ^ acting on galaxies as 
function of clu ster-ce ntric radius and orbital velocity 
(Haines et al. l201ld ). For the cluster-centric radius 
we assume the projected distance of 0.995 Mpc. As 
a lower limit for the orbital velocity, we assume the 
galaxy's line-of-sight velocity relative to the cluster's 
systemic velocity V\ os ~830 km s , For the upper 
limit we assume the escape velocity at the cluster- 
centric distance which, from the Henrquist model for 
A 3558, is ~2600kms~ 1 . At 0.995 Mpc from the clus- 
ter centre picm — 1.3T0 -28 gem -3 and the ram pres- 



sure ranges from ~9- 10 dyncm to ~9T0 



-12 



This is the regime of weak to mode rate r am pres 
sure as defined by Roediger & Hensler (j2005h who, by 
means of high resolution 2D hydrodynamical simula- 
tions, demonstrated that ram pressure effects can be 
observed over a wide range of ICM conditions. In par- 
ticular, in high density environments RPS severely 
truncates the gas disc of L* galaxies, while in low 
density environments, where moderate ram pressure 
is for eseen, their gas disc is cl early disturbed and 
bent (Roediger & Henslerl . I2OO5I ) as we observe. The 
gas discs of these galaxies can be truncated to 15- 
20kpc in the first 20-200 Myr of RPS. To check if 
the expected ram pressure is able to strip the gas 
from the disc of SOS 114372, we estimate the gravi- 
tational restoring force per unit area (d(j)/dz)cr gas = 
2ttG 'o disco 'gas, with a s tar and a disc the star and gas 
surface density, respectively. From the -fT-band pho- 



tain M/Lk ~2, therefore the stellar mass surface 
density is 8.8M pc -2 . Assuming a typical a gas ~ 
IOMqpc -2 , the restoring force per unit area is 1.7 • 
10 _12 dyncm~ 2 which compared with the estimated 
ram pressure (9-10 _13 -9T0 -12 dyncm -2 ) shows that 
it can be effective in removing the gas from the disk 
if Vorb > 1100 km s . We assume that the galaxy is 
coming in radially, such that the projected cluster- 
centric distance is the same as the actual 3D one. If 
the galaxy is at a larger 3D distance, then the lo- 
cal density will be lower, as will the wind velocity, 
and thus the ram pressure could be overestimated al- 
though a moving ICM and/or local overdensities can 
enhance it. 

Crowl & Kenney (2006. 120081 ) conducted a sur- 
vey of ten Virgo galaxies with the SparsePak IFS. 
They found clear evidence of RPS-induced star for- 
mation within the truncation radius and a passive 
population beyond it. Their estimates of the strip- 
ping timescales showed that Virgo galaxies can be 
stripped both within and beyond the cluster core. 
This is direct evidence that the ICM is neither static 
nor homogeneous. That a galaxy's ISM can be af- 
1 felted by the environment also in intermediate- to 
^ iow-density regions (0.6-1.0 Mpc from the cluster cD 
galaxy) in the Virg o cluster was also concluded by 



Chung et al.l (|2007l ) observing the HI tails of seven 



spiral galaxies. 

A non homogeneous ICM is expected for non viri- 
alized mergi ng and post-merging clusters. Dynami- 
cal analyses (jBardelli at al. . 19981 ) and X-r ay obs er- 
vations (Bardelli et al. [l996j; Akimoto et al. 120031 ) of 
A 3558 indicate this cluster is the result of a cluster- 
cluster merging and/or a group infall. More recently, 



combi ning XMM-Newton and Chandra data. iRossetti et al 



(2007|) were able to detect a cold front in the A 3558 
ICM probably caused by the sloshing of the core in- 
duced by the perturbation of the gravitational poten- 
tial associated with a past merger. This cold front is 
expanding, perturbing the density and the tempera- 
ture of the ICM. From their Fig. 3 it is clear that the 



cold front is moving in the direction of SOS 114372 
tometry, we derive the inclination- and dust attenuation- r • / nr,m ■ ■ x- ^ 

and that it is ~5arcmm (~ 300 kpc m projection) 

corrected surface brightness at 20 kpc (where the gas 
disc is bent) of 4.4L Q pc~ 2 . From the i^-band lu- 
minosity and the stellar mass (see Table [1]) , we ob- 



from the galaxy. Unfortunately, the X-ray observa- 
tions do not include the region of SOS 114372, which 
is planned to be observed in future surveys (see Hu- 
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to RPS despite its large distance from the cluster cen- 
tre. 
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N-body/hydrodynamical simula- 
tions of ram-pressure stripping 



daveri, Bozkurt & Ercan 2O10h . rotation of the galaxy breaks one important symme- 

Taking into account all the above arguments, we try that lies at the ground of the applicability of that 
conclude that the environment of our galaxy is favourabkcheme. The second reason is that not only the di- 
rection of the outflow does not indicate the direction 
of the wind, but also the direction of the wind with 
respect to the disc cannot be assumed a priori (see 
Appendix A). We therefore considered all the possi- 
ble values for the angle /3 between the galaxy rotation 
axis and the ICM wind. We ran simulations from 
RPS close to edge-on (ft = 85°) to close to face-on 
(/3 = 15°) at intervals of 10° (there was no need to ex- 
plore values closer to face-on or to edge-on) . For what 
concerns the ICM wind velocity, we set as the lower 
limit the galaxy's line-of-sight velocity relative to the 
cluster's systemic velocity, V\ os ~830 km s ^ , and as 
the upper limit the escape velocity, ~2600kms~ 1 
(derived from the Hernquist model for A 3558). 

We explored the following values for the ICM 
wind velocity (V wind =830, 1100, 1400, 1700, 2000, 
2200, 2400, 2600 kms" 1 ). We also explored two dif- 
ferent disc scale heights corresponding to 10% and 
20% of the radius r<j. In summary, our grid of mod- 
els includes 8 values for /?, 8 values for V w i n d arid 2 
values for the disc scale height. Our approach was to 
run low-resolution simulations for the whole range of 
parameters and then select the most promising cases 
for the high-resolution simulations. 

The simulations were done with the cosmological 
N-body/hydrodyna mic code GADG ET- 2 developed 
by Volker Springel (jSpringell . 120051') . The hydrody- 



To confirm the scenario of RPS and to constrain its 
physics, we performed N-body/hydrodynamical sim- 
ulations. Our aim was to reproduce as closely as 
possible the velocity field of the gas. We remark that 
in the real system there are factors that cannot be 
accounted for in a simulation, such as the spatial dis- 
tribution of the dust, which may affect the observed 
kinematics and distribution of the gas. The interac- 
tion of the ICM with the ISM may also be compli- 
cated by their dumpiness, not to mention the action 
of the bar. We therefore don't expect the simula- 
tions to provide a description of the details of the gas 
kinematics, but aim to obtain an overall agreement 
with the observations helping us to constrain the ICM 
wind angle and velocity and the 'age' of RPS. 

9.1 Simulations 

An important aspect of the use of simulations is that, 
at least in the early phases of RPS, the direction of 
the ICM wind cannot be inferred from the apparent 
direction of the gas outflow, because the motion of 
the gas is a complex combination of the effects of the 
ram pressure and the gravitational fie ld of the host 
galaxy (see a lso Roedigger & Briiggen 
et al. ( 19981 ) introduced a geometrical scheme (see 
their Fig. 1) to study the jets in Seyfert galaxies. 
In that scheme, the angle between the jet and the 
normal to the disc is constrained by the apparent 
jet outflow angle and the inclination of the galaxy. 
In a different framewo rk, t hat scheme was adopted 
by Abramson et al. ( 201ll ) to study the geometry 
of RPS, substituting the jet direction with the ICM 
wind direction. 

Unfortunately, we cannot adopt Clarke et al.'s 
scheme for two reasons. The first reason is that the 



namic part follow s an SPH scheme (|Gingold Monaghan . 
19771 : iLucvl . Il977h while the gravitational interaction 



is calculated using a tree code for the short r ange 
20061 ) Clarke f° rce as nrs * introduced by Barnes Hut ( 19861 ) and 
treePM code based on Fourier techniques for the long 
range forces. Furthermore, o ur version of GADG ET- 

' 19961 ) and 



2 includes radiative cooling (jKatz et al 



recipes for star formati on, stellar feedback and gala c- 
tic winds according to Springel Sz Hernquist ( 20031 ). 

To simulate galaxies undergoing RPS, we u s e the 

same si mula tion setup as iKapferer et all (j2009bh . iKronberger et al 
(|2008ah and lSteinhauser et al.1 (|2012h . First, we gen- 
erate a model galaxy matching the observed prop- 
erties of SOS 114372 (see Sect. [2]). The initial condi- 
tions are calculated ac cording to 



based on the work of Mo et al 



Springel et al.1 (|2005h 



,P 
(1 



1998). Our model 
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Table 2: The model galaxy. Properties of the initial conditions of the model galaxy used. The particle 
numbers, mass resolution and cumulative mass are shown for disc scale length corresponding to 10% of the 
disc radius both for low- and high-resolution simulations (LR, HR). 





# of particles 


mass resolution 


total mass 






[M & / particle] 


[M Q ] 


HR 






DM halo 


3 x 10 5 


3.7 x 10 6 


1.1 x 10 12 


gas disc 


2 x 10 5 


2.9 x 10 4 


5.9 x 10 9 


stellar disc 


2 x 10 5 


2.6 x 10 5 


5.3 x 10 10 


stellar bulge 


1 x 10 5 


1.7 x 10 5 


1.7 x 10 10 


LR 






DM halo 


3 x 10 4 


3.7 x 10 7 


1.1 x 10 12 


gas disc 


2 x 10 4 


2.9 x 10 5 


5.9 x 10 9 


stellar disc 


2 x 10 4 


2.6 x 10 6 


5.3 x 10 10 


stellar bulge 


1 x 10 4 


1.7 x 10 6 


1.7 x 10 10 



galaxy includes a dark matter halo, a stellar disc, 
a gaseous disc and a stellar bulge. The disc scale 
length of the galaxy is directly related to the angu- 
lar momentum of the disc with the spin as a free 
parameter. The stellar and gaseous disc have expo- 
nential surface density profile s. For the s tellar bulge 
and the dark matter halo a iHernauistl (|l990h pro- 
file is adopted. Important parameters and proper- 
ties adopted for our model galaxy are shown in Ta- 
ble [2] for both the low- and high-resolution simula- 
tions. The total mass is .Mtot=1.2-10 M© and the 
total (dynamical) and baryonic mass within 20kpc 
are M dyn =2.1-10 11 M & (see Sect.E2D and A4*=7.0-10 10 
respectively. Since the initial gas mass of SOS 114372 
is not known, we assume that the gaseous phase con- 
tains 10% of the mass of the disc. The initial disc 
scale height was set to either 10% or 20% of the 
disc scale length. The disc-to-bulge mass ratio of 
0.32 adopted for the model galaxy was derived from 
the measured B/D flux ratio of 0.44 taking into ac- 
count that ~30% of the light from the bulge is due 
to very young stars contributing only 1% to its mass 
(Sect. [72}. 

Since the model galaxy needs some time to settle 
to avoid numerical artefacts, it is evolved in isolation 



for 1 Gyr before the ICM wind is applied. Further- 
more, the parameters for generating the model galaxy 
were chosen carefully to obtain a stable galaxy. As 
also found by lMoster et al.1 (1201 ll ). the resolution plays 
an important role. As we found in the low-resolution 
runs in which the galaxy consists of a significantly 
smaller number of particles than the nominal 7.5 x 10 5 
value, both the stellar and the gaseous disc were 
puffed up, leading even to warps although the galaxy 
is evolving in isolation. With the chosen parameters 
the effect of the resolution was less dramatic and neg- 
ligible on the morphology of the stripped gas, used 
Mg^he selection of the models which better reproduce 
the observations. 

To simulate the prevailing environment of 
SOS 114372, the model galaxy is then exposed to an 
external wind, simulating the ICM. At variance with 
many other RPS simulations, the model galaxy itself 
is moving and the ICM is stationary. The gas par- 
ticles representing the ICM fill a simulation domain 
with periodic boundary conditions in order to keep 
the ICM density stable during the simulation. The 
mass resolution of the ICM particles is the same as 
for the ISM particles to avoid numerical artefacts. 
As we use a box size of 850 kpc, a large number of 
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particles would be needed leading to a very long com- 
putation time. For this reason, only the inner part 
of the cube with size 200 x 200 x 850 kpc 3 , where the 
galaxy flies through, has the same mass resolution as 
the ISM. On the outside, particles with a hundred 
times higher mass are used. For an ICM density of 
Picv\ = 1.3 x 10 -28 gcm -3 ( Sanderson &: Ponmanl . 



2010). 3 .1 million gas particles are needed. 



The models were run from the onset of the ram 
pressure to 100 Myr after. During this time the ram 
pressure is kept constant and its onset is immediate. 

9.2 Comparison with the observations 

To compare the models with the data we first ac- 
counted for the constraints on the viewing directions. 
The first constraint is given by the inclination an- 
gle i=82° of the galaxy. The other constraint comes 
from the requirement that the spatial velocity of the 
galaxy V ga i projects into the observed Vi os . This 
translates into the requirement that the angle (ft be- 
tween the motion of the galaxy and the line of sight is 
given by cos (ft = Vi os /V ga i. As it is shown in the Ap- 
pendix, the consequence of these constraints is that 
some models can be always observed from four direc- 
tions, others from two, and others from two but only 
for high values of V ga i. The issue of the viewing di- 
rections is treated in detail in the Appendix where we 
also describe the process leading from the 3-D models 
provided by the simulations to the 2-D fields compa- 
rable with our data. In brief, to compare the models 
with the observations we account for the seeing ef- 
fect associating to each particle of the simulation a 
Gaussian distribution with the FWHM equal to the 
seeing, and in each surface element of the plane of 
the sky we compute the surface density of particles 
£ as the sum of the contributions from the Gaussians 
and the radial velocity as the weighted average of the 
individual velocities. In the following, a case corre- 
sponding to a particular projection of a simulation 
(characterized by wind velocity, wind angle and time 
past since the onset of RPS) will be referred to as 
'case' or 'model', without distinction. 

For what concerns the spatial distribution of the 
gas, we should keep in mind that while the simula- 
tions show the velocity field of the whole gas, our 



data refer to the ionized gas. However, it is reason- 
able to assume that the neutral and ionized gas are 
mixed. In their study of the ram-pre ssure stripped 
galaxy NGC4330, Abramson et al. (|201ll ) showed 
that the Ha flux is truncated at the same radius on 
the SE side as the HI emission which extends slightly 
furthe r in the trailing side, and Tonnesen & Bryan 
(|2O10h found that in the tail of the Ha emission fol- 
lows the neutral gas. 

The first selection was made on low-resolution 
models. To start the selection, we compared the spa- 
tial distributions and the velocity fields of the models 
with the observations as explained below. 

The Ha surface brightness was first normal- 
ized to 5xl0 _18 ergs— lcm~ 2 A— larcsec -2 so that at 
our 5a detection limit F=l. Then, for each case, we 
computed the quantity x 2 = J2{ F ij~^ij) 2 / (N pix -T), 
where £ is the particle surface density of the model 
and Npi x is the number of pixels where x 2 is evalu- 
ated. The sum was computed, in a first iteration, over 
the intersection between the observed area and the 
area where £ >0. In this computation we masked out 
the two regions with very high star- formation rates, 
i.e. the nuclear region and the region 12 kpc SW from 
the nucleus. We then ranked the cases in increasing 
order of x 2 an d by visual inspection verified that the 
ranking in x 2 reflects faithfully the consistency of the 
model gas distribution with the data. 

The observed gas is confined in an area limited 
by our detection limit of F=l, while the £ of the 
models may assume arbitrarily small values. To set 
a lower limit to £ consistent with the observational 
limits, we derived, for each case, the mean a of the 
ratios Sjj/i^,- computed over the flux range .F=l-10 
(i.e. from our flux limit to 10 x higher), and derived a 
mean value of a averaging over the first ranked cases. 
We found a=1.21±0.23, with the low dispersion re- 
flecting the stability of a over ~ 100 different models. 
By definition, a is the particle surface density corre- 
sponding to our lowest detection value, i.e. to the 
"borders" of our data, and as such it was adopted as 
the cutoff value to compare the LR models with the 
observations. For the HR simulations, having 10 time 
higher particle density, we used a cutoff 10 x higher 
than this (a'=12.1) . 

We then run a second iteration of the computa- 
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Figure 18: Example of a grid of low-resolution models. The wind angle /3 (in degrees) changes in the vertical 
direction and the time since the onset of RPS (in Myr) in the horizontal direction as indicated. The selected 
models are marked with a light circle in the top-right of the thumbnails. Notice that all of the models are 
shown for brevity only in one projection(case B, see Appendix). V w i n d is fixed at 1400 kms" 1 . The colour 
bar gives the values of radial velocity in kms -1 . Notice the progressive bending of the disc with time and 
the progressive truncation of the upper-left disc with /3. 
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Figure 19: Analogous to Fig. \18\ but with the wind velocity V w i n d in place of the time. The time is fixed at 
60 Myr. 
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Figure 20: Left: One of the best models ( j3 = 45°, 14,^=1400 km s _1 , t=60Myr; left) is compared to the 
observed radial velocity field. The arrow shows the direction of the ICM wind. The velocity scale is to the 
right. Right: The contours of Ha flux measured with WiFeS are superimposed to the gas density of the 
same model. 



tion of x 2 ) where we set S > a to cut the models 
consistently with the data. We also considered the 
size of the intersection between the observed area 
and the models, and accepted only the models which 
cover at least 80% of the observed area. In this way 
we could reject models whose spatial gas distribu- 
tion is significantly different from the observed one. 
These are the cases when a substantial fraction of the 
gas is stripped from the disc, like in advanced phases 
of RPS (t>70 Myr) or in early phases (t<20 Myr), 
where the gas in the models is still confined in the 
disc. Only 38% cases survived this selection. 

A similar ranking was applied to the velocity fields, 
defining X 2 = E(K? " V°A 2 /(N pix - 1) 



where VJ? 



and V°j are the model and observed radial velocity 
fields respectively. In this case the ranking accord- 
ing to increasing values of x 2 is more prone to 'mis- 
takes' than before, and actually the preferred cases 
are those at the extremes of time: t>80 Myr or t<20 
Myr. The reason is that those models spatially over- 
lap with the data only in the inner regions, where 
the velocity field of the gas is less perturbed. To 
remove the contamination from these cases, we had 
to exclude models with overlap less than 90%, which 
left only 5% of them. Also in the case of velocity 
fields we verified by visual inspection that the rank- 
ing reflects the consistency between data and models. 
The cases that survived this selection were further in- 
spected visually to look for manifest inconsistencies 



between models and data. 

Figures [18] and [19] show some models in two pa- 
rameter grids (/3,t and (3,V w i n d respectively). Some 
selected models are marked with a circle on the upper- 
left corner of the thumbnails. From the low-resolution 
simulations, we selected the models with /3=45° and 
(3=55° with ^^=1400 km s" 1 and 1700 km s" 1 , and 
disc scale height 10% of the scale length, and carried 
out high-resolution simulations with these parame- 
ters. Lower and higher velocity values at the same f3 
values are not able to produce the observed gas out- 
flowing or truncate the gas disc at smaller radii than 
observed, respectively (see Fig. [T9]) . At a certain ve- 
locity value, (3 values closer to the edge-on RPS tends 
to truncate the NE galaxy side sharply respect to the 
observations, while (5 values closer to the face-on RPS 
presents blueshifted gas in the NW side of the disc in 
contrast with what observed (see Figs. [181 and [7|) . 

The comparison of the observed velocity field with 
the high-resolution simulations leads to a further re- 
striction of the parameter space. The best final model 
is shown in Fig. [20] (left panel) together with the ob- 
served velocity field. The projection on the plane of 
the sky of the ICM wind is also shown. This model 
corresponds to (3 = 45°, 14,^=1400 km s _1 and t = 
60 Myr. Figure [20] (right panel) compares the gas 
density of the same model to the observed Ha flux 
distribution. The most external isophote corresponds 
to our 5a detection limit (F = 1) while the density 
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of the model is consistently cut at £ = 12.1 parti- 
cles per arcsec 2 (see above). We notice a remarkable 
agreement between the model and WiFeS observa- 
tions. In the southwestern part of the disc not cov- 
ered by WiFeS data the model is less constrained 
by the Ha imaging (cf. Fig. [6] right panel) which is 
slightly shallower than the IFS observations. 

Our data are therefore consistent with a wind 
angle j3 ~ 45°, such that the RPS is taking place 
at an angle intermediate between face-on and edge- 
on. The wind velocity is 1400 km s -1 (ram pressure 
~ 9 • 10 12 dyncm -2 ). The only orientation consis- 
tent with the observations corresponds to the case in 
which we are observing the galaxy disc from the side 
of impact of the ICM wind, in such a way that the 
SE side of the disc is the furthest from us. All se- 
lected models correspond to the disc height scale of 
0.1 xr^. The range of the parameters is very narrow. 
Models with f3 < 45° or /3 > 55° are clearly less con- 
sistent with the observations, as are the models with 
Vwind > 1700 km s -1 or V wind < 1400 km s -1 . For 
what concerns the time since the onset of the RPS, 
the best models are those at t~60 Myr. We notice 
that the time scale of the RPS suggested by the sim- 
ulations is indicative, but the observational evidence 
of the outflowing gas following the rotation of the 
galaxy shows that the stripping process has started 
recently. Considering the adopted value of the ICM 
density and the estimated wind velocity we obtain a 
ram pressure of 2.5-10~ 12 dyncm -2 . 

10 Summary and Conclusions 

We have presented IFS observations of the bright spi- 
ral galaxy SOS 114372, located at 1 Mpc from the 
centre of the rich cluster A 3558 in the Shapley su- 
percluster core. IFS data demonstrate the presence 
of a one-sided ionized gas stream extending 13kpc 
in projection from the galaxy disc. Complementary 
narrow-band Ha imaging shows multiple compact knots 
of Ha emission often connected to faint filamentary 
strands which extend back to the galactic disc, while 
the -fT-band image shows a symmetric disc. These 
data, complemented by UV, FIR and deep optical- 
NIR (B, R, K) imaging, have allowed us to inves- 



tigate the gas and stellar kinematics as well as the 
properties of the gas and stellar populations which 
all together indicate that this galaxy is being affected 
by ram-pressure stripping. 

The observations have been interpreted by means 
of emission-line diagnostics, stellar population syn- 
thesis and photoionization models and N-body/hydrodynamical 
simulations. The main results of our analysis are the 
following. 

i) The velocity field of the gas is complex and, al- 
though still dominated by the potential of the 
galaxy, is affected by significant perturbations 
even in the inner parts of the disc. The ve- 
locity dispersion field is characterized by max- 
ima (^ 100 km s _1 ) localized in the periphery of 
the disc. The velocity dispersion maintains rel- 
atively high values (;> 40kms~ 1 ) all over the 
field. At the same time, the stellar velocity 
field appears to be symmetric with respect to 
the centre, with the typical shape of a barred 
galaxy. 

ii) All the gas in the galaxy is contaminated by some 
fraction of shock excitation (0.05 / 0.1) as 
shown by the shock and photoionization mod- 
els. The galaxy regions associated to the out- 
flowing gas and three regions on the opposite 
side of the disc are dominated by shock-excited 
gas (0.4 ^ / ^ 0.8), and all but one are strongly 
attenuated by dust. In these regions we also 
measure the highest gas velocity dispersions. 

Hi) The current star formation rate, measured through 
the Ha luminosity, from IFS observations, amounts 
to 7.2±2.2 Moyr -1 and shows two prominent 
peaks in the centre and at 12 kpc to the SW, 
which contribute 30% and 20% to the total SFR, 
respectively. This picture is fully consistent 
with that measured from the 24^m flux. The 
analysis of the stellar populations led to the 
identification a starburst characterized by a ~5 x 
increase in the SFR over the last ~100Myr in 
the SW disc, a short starburst (~ 300 Myr) just 
quenched in the NE disc and in the galaxy cen- 
tre a predominant population of old (>lGyr) 
stars together with newly (<^100 Myr) formed 
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stars. The intense nuclear star formation (~ 2 
Mq yr _1 ) is most probably associated with the 
presence of the bar, whose effect is to concen- 
trate the gas to wards the cent r e thu s triggering 
star formation ( Ellison et all 2011 ; Oh et al. 
2012L and references therein). 



most luminous (Mg ~-21.5) and distant (2,43558 =0-048) 
galaxy showing ongoing RPS studied in detail (c.f. 
galaxies in Virgo, Coma and A 1367 clusters; e.g. 
Chun g et a l. 120091 and refe rence therein; Abramson 



iv) Gas flowing out of the discs of galaxies is ob- 
served also in case of tidal interactions or galac- 
tic winds. Among the closest galaxies to 
SOS 114372 we do not find any likely compan- 
ion responsible for the gas outflow. In fact, 
SOS 112392 and SOS 114493 present regular mor- 
phologies with no signs of interaction. The other 
candidate for the tidal interaction, the faint 
galaxy SOS 115228 with unknown redshift, shows 
a starburst which could be related to an inter- 
action with SOS 114372, but its very low mass 
as well as the relative positions of these two 
galaxies cannot justify the observed gas out- 
flow. We also considered the possibility of a 
galactic wind, but the morphology of the out- 
flow does not agree with the bipolar outflows 
expected in that case. On the other hand, the 
high velocity dispersion and the gas moving to- 
wards the observer in a triangular region along 
the minor axis SE from the nucleus is suggestive 
of the semi-cone of a galactic wind produced by 
the nuclear starburst with the other semi-cone 
obscured by dust. In this case the wind would 
be also affected by the RPS and dispersed on a 
short time scale. 

The RPS scenario for SOS 114372 is also sup- 
ported by ad hoc iV-body/hydrodynamical simula- 
tions. They allowed us to constrain the angle between 
the rotation axis and the ICM wind, f3 ~ 45°, the 
wind velocity, V w i n d ~ 1400kms _1 , and the "age" 
of the gas stripping, t~60Myr. The fact that the 
ram-pressure stripping started 60Myr ago is consis- 
tent with the burst of star formation in the SW disc 
starting less than 100 Myr ago according to the stel- 
lar population modeling. We are thus led to ascribe 
this starburst to the compression of the ISM caused 
by the ram pressure. 

In this work we have analyzed the case of an L* 
galaxy affected by ram-pressure stripping. It is the 
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SOS 114372 is located far away from the cluster cen- 
tre, in an environment characterised by a relatively 
low density with respect to the cluster's core. This 
adds a piece of evidence to the fact that RPS is acting 
more efficiently than previously foreseen and also out- 
side of th e cluster cores as als o observed in the Virgo 
cluster by Chung et al. ( 20071 ) . This new understand- 
ing of the RPS candidates supports the view that this 
mechanism is the pri ncipal transformation proces s to 
turn spirals into SOs (jKormendv &: Bender . 2012 ) al- 
though helped by other processes affecting the struc- 
ture of the galaxies. 

We plan to analyse the supercluster sample ob- 
served with WiFeS in forthcoming papers providing 
other examples of the influence of the environment 
on galaxy evolution. 

Acknowledgments 

This work was carried out in the framework of the col- 
laboration of the FP7-PEOPLE-IRSES-2008 project 
ACCESS. CPH acknowledges financial support from 
STFC. RJS was supported under STFC rolling grand 
PP/C501568/1 "Extragalactic Astronomy and Cos- 
mology at Durham 2008-2013" . Dopita acknowledges 
the support from the Australian Department of Sci- 
ence and Education (DEST) Systemic Infrastructure 
Initiative grant and from an Australian Research Coun- 
cil (ARC) Large Equipment Infrastructure Fund (LIEF) 
grant LE0775546 which together made possible the 
construction of the WiFeS instrument. Dopita would 
also like to thank the Australian Research Council 
(ARC) for support under Discovery project DP0984657. 
DS and SS acknowledge the Unilnfrastrukturprogramm 
des BMWF Forschungsprojekt Konsortium Hochleis- 
tungsrechnen and the Austrian Science Foundation 
FWF through grant W1227 (FWF DK-CIM) and 
thanks Volker Springel for providing the simulation 
code GADGET-2 and A. Bohm for profitable discus- 
sions. The authors wish to thank the staff at the Aus- 



38 



tralian National University 2.3m telescope for their 
assistance for WiFeS observations. PM, GB and CPH 
thank the ANU-RSAA for hospitality. PM and GB 
thank C. J. Farage for helpful discussions about the 
use of WiFes pipeline. CPH thanks Sylvain Veilleux 
for his assistance in the acquisition and calibration 
of the Magellan/MMTF Ha imaging and the techni- 
cal staff at the Las Campanas Observatory for their 
support during the observations. A special thanks 
to A. J. R. Sanderson for providing the X-ray pro- 
file of A 3558 and to K.A. Pimblett and A. W. Gra- 
ham for their comments which helped to improve the 
manuscript. Figures 1 and 5 were prepared with the 
IDL code provided by Lupton et al. (2004) and avail- 
able at http:/ /cosmo. nyu.edu/hogg/visualization/rgb/ 
The authors thank the referee B. Vollmer for his con- 
structive comments and suggestions. 

References 

Abadi M. G., Moore B., & Bower R. G., 1999, MN- 
RAS, 308, 947 

Abramson A., Kenney J. D. P., Crowl H. H. et al., 
2011, AJ, 141, 164 

Akimoto F., Kondou K., Furuzawa A., Tawara Y., 
Yamashita K., 2003, ApJ, 596, 170 

Allen D. P., Norris R. P., Meadows V. S., Roche P. F., 
1991, MNRAS, 248, 528 

Ascasibar Y. & Diego J. M., 2008, MNRAS, 383, 369 

Athanassoula E. & Misiriotis A., 2002, MNRAS, 330, 
35 

Baes M., Davies J. I., Dejonghe H., et al., 2003, MN- 
RAS, 343, 1081 

Baldry I. K., Balogh M. L., Bower R. G., et al., 2006, 
MNRAS, 373, 469 

Balogh M. L., Navarro J. F., & Morris S. L., 2000, 
ApJ, 540, 113 

Bamford S. P., Nichol R. C, & Baldry I. K., 2009, 
MNRAS, 393, 1324 



Bardelli S., Zucca E., Malizia A., et al., 1996, A&A, 
305, 435 

Bardelli S., Pisani A, Ramella M., Zucca E., 
Zamorani G., 1998, MNRAS, 300, 589 

Barnes J. & Hut P., 1986, Nature, 324, 446 

Barnes J. E. & Hernquist L. E., 1991, ApJ, 370, 65L 

Bekki K., 2001, ApJ, 546, 189 

Bekki K., Couch W. J., & Shioya Y., 2002, ApJ, 577, 
651 

Bekki K., 2009, MNRAS, 399, 2221 

Boselli A. & Gavazzi C, 2006, PASP, 118, 517 

Bosma A., 1981, AJ, 86, 12 

Blanton M. R., Eisenstein D., Hogg D. W., et al., 
2005, ApJ, 629, 143 

Bureau M. & Athanassoula L., 2005, ApJ, 626, 159. 

Butcher H. & Oemler Jr. A., 1984, ApJ, 284, 426 

Byrd G. & Valtonen M., 1990, ApJ, 350, 89 

Cardelli J. A., Clayton G. C, & Mathis J. S., 1989, 
ApJ, 345, 245 

Cappellari M. & Copin Y., 2003, MNRAS, 342, 345 

Chemin L., Cayatte V., Balkowski C, et el., 2005, 
A&A, 436, 469 

Chung A. & Bureau M., 2004, AJ, 127,3192 

Chung A., van Gorkom J. H., Kenney J. D. P., 
Vollmer B., 2007, ApJ, 659L, 115 

Chung A., van Gorkom J. H., Kenney J. D. P., Crowl 
H., Vollmer B., 2009, AJ, 138, 1741 

Cid Fernandes R., Mateus A., Sodre, Stasihska G., 
Gomes J. M., 2005, MNRAS, 358 363 

Conselice C. J., Bundy K., Ellis R. S, Brichmann J., 
2005, AJ, 628, 160 



Cortes J. R., Kenney J. D. P., & Hardy E., 2006, AJ, 
131, 747 



39 



Cortese L., Marcillac D., Richard J., et al., 2007, 
MNRAS, 376, 157 

Cortese L., Davies J. I. Pohlen M., Baes M., et al., 
2010, A&A, 518, 49 

Courteau S., Andersen D.R., Bershady M.A., 
MacArthur L.A., Rx H.-W., 2003, ApJ, 594, 208 

Clarke C. J., Kinney A. L., & Pringle J. E., 1998, 
ApJ, 495, 189 

Cowie L. L. & Songalia A., 1977, Nature, 266, 501 

Crowl H. H., Kenney J. D. P., van Gorkom J. H., 
Vollmer, B., 2005, A J, 130, 65 

Crowl H. H. & Kenney J. D. P. 2006, ApJ, 649L 75 

Crowl H. H. & Kenney J. D. P. 2008, AJ, 136, 1623 

Dale D. A., Giovaneli R., Haynes M. P., Hardy E., 
Campusano J. E., 1999, ApJ, 118, 1468 

Dale D. A., Giovaneli R., Haynes M. P., Hardy E., 
Campusano J. E., 2001, AJ, 121, 1886 

de Naray R. K., Zagursky M. J., McGaugh S. S., 
2009, AJ, 138, 1082 

Diehl S. & Statler T. S., 2006, MNRAS, 368, 497 

Dopita M. A. & Sutherland R. S., 2003, Astrophysics 
of the Diffuse Medium, Springer- Verlag Berlin Hei- 
delberg 2003 

Dopita M. A., Hart J., McGregor P., Oates P., Blox- 
ham, G., Jones D., 2007, Ap&SS, 310, 255 

Dopita M. A., Rhee J., Farage C. et al., 2010, Ap&SS, 
327, 245 

Dressier A., 1980, ApJ, 236, 351 

Dressier A., Oemler A. Jr., Couch W. J., et al. 1997, 
ApJ, 490, 577 

Dressier A., Smail I., Poggianti B. M., et al. 1999, 
ApJS, 122, 51 

Dressier A., Oemler A. Jr., Poggianti B. M., et al. 
2004, ApJ, 617, 867 



Driver S.P., Popescu C.C., Tuffs R.J., et al., 2008, 
ApJ, 678, L101 

Ellison S.L., Preethi N., Patton D.R., et al, 2011, 
MNRAS, 416, 2182 

Emsellem E., Cappellari M., Krajnovic D., et al., 
2007, MNRAS, 379, 401 

Farage C. L., McGregor P. J., Dopita M. A., Bicknell 
G. V., 2010, ApJ 724, 267 

Fathi K., van de Ven G., Peletier R.F., et al., 2005, 
MNRAS, 364, 773 

Fischera J. k Dopita M. A., 2005, ApJ, 619, 340 

Fossati M., Gavazzi G., Boselli, Fumagalli M., 2012, 
ASA, 544, 128 

Fujita Y. & Nagashina M., 1999, ApJ, 516, 619 

Gavazzi G., Boselli A., Mayer L., et al., 2001, 563L, 
23 

Gingold R. A. & Monaghan J. J., 1977, MNRAS, 181, 
375 

Giovanelli R. & Haynes M. P., 2002, ApJL, 571, 107 

Graham A. W. & Worley C. C, 2008, MNRAS, 388, 
1708 

Grevesse N., Asplund M., Sauval A. J., Scott P., 
2010, Ap&SS, 328, 179 

Gunn J. E. & Gott J. R. Ill, 1972, ApJ, 176, 1 

Haines C. P., Merluzzi P., Mercurio A., et al. 2006, 
MNRAS, 371, 55 

Haines C. P., Gargiulo A., La Barbera F., Mercurio 
A., Merluzzi P., Busarello G., 2007, MNRAS, 381, 
7 

Haines C. P., Smith G. P., Egami E., et al., 2009, 
ApJ, 704, 126 

Haines C. P., Busarello G., Merluzzi P., et al., 2011a, 
MNRAS, 412, 127, Paper II 

Haines C. P., Busarello G., Merluzzi P., et al., 2011b, 
MNRAS, 412, 145, Paper III 



40 



Haines CP., Merluzzi P., Busarello G., et al., 2011c, 
MNRAS, 417, 2831 

Haynes M. P., Giovanelli R., & Kent B. R., 2007, 
ApJ, 665L, 19 

Hernquist L., 1990, ApJ, 356, 359 

Hoopes C.G., Heckman T.M., Strickland D.K., et al., 
2005, Apj, 619, L99 

Hundaverdi M., Bozkurt M., & Ercan N., 2010, ASP 
Conference Series, Vol. 424, 330 

Irwin J. A. & Sarazin C. L., 1996, ApJ, 471, 683 

Jachym P., Koppen J., Palous J., Combes F., 2009, 
A&A, 500, 693 

James P. A., Prescott M., & Baldry I. K., 2008, A&A, 
484, 703 

Jimenez- Vicente J., Mediavilla E., Castillo-Morales 
A., Battener E., 2010, MNRAS, 406, 181 

Kapferer W., Kronberger T., Ferrari C, Riser T., 
Schindler S., 2009a, MNRAS, 389, 1405 

Kapferer W., Sluka C, Schindler S., Ferrari C, 
Ziegler B., 2009b, A&A, 499, 87 

Katz N., Weinberg D. H, Hernquist L., 1996, ApJS, 
105, 19 

Kauffmann G., Heckman T. M., Tremonti C, et al. 

2003, MNRAS, 346, 1055 

Kauffmann G., White S. D. M., Heckman T. M., et 
al., 2004, MNRAS, 353, 713 

Kenney J. D. & Koopmann R., 1999, AJ, 117, 181 

Kenney J. D. P., van Gorkom J. H., & Vollmer B., 

2004, AJ, 127, 3361 

Kenney J. D. P., Tal T., Crowl H. H., Feldmeier J., 
Jacoby G. H., 2008, ApJ, 687L, 69 

Kennicutt R. C. Jr., 1983, ApJ, 272, 54 

Kennicutt R. C. Jr., Tamblyn P., Congdon C. E., 
1994, ApJ, 435, 22 



Kennicutt R. C. Jr., 1998, ApJ, 498, 541 

Kewley, L. J., Heisler, C. A., Dopita, M. A., Lums- 
den, S. 2001, ApS, 132, 37 

Kewley, L. J., Groves, B., Kauffmann G., Heckman 
T., 2006, MNRAS, 372, 961 

Kim D.-W., Kim E., Fabbiano G., Trinchieri G., 
2008, ApJ, 688, 931 

Kimura H., Mann I., & Jessberger E. K., 2003, ApJ, 
582, 846 

Koleva M., Prugniel Ph., Ocvirk P., Le Borgne D., 
Soubrian C, 2008, MNRAS, 385, 1998 

Koleva M., Prugniel Ph., Bouchard A., Wu Y., 2009, 
A&A, 501, 1269 

Koopmann R. A., Giovanelli R., Haynes M. P., et al., 
2008, ApJ 682L, 85 

Kormendy J. & Bender R., 2012, ApJSS, 198, 2 

Kregel M., van der Kruit P. C, & Freeman K. C, 
2004, MNRAS, 351, 1247 

Kronberger T., Kapferer W., Schindler S., et al., 
2006, A&A, 458, 69 

Kronberger T., Kapferer W., Schindler S., Ziegler 
B. L., 2007, A&A, 473, 761 

Kronberger T., Kapferer W., Ferrari C, Unter- 
guggenberger S., Schindler S., 2008a, A&A, 481, 
337 

Kronberger T., Kapferer W., Unterguggenberger S., 
Schindler S., Ziegler B. L., 2008b, A&A, 483, 783 

Kroupa P., 2001, MNRAS, 322, 231 

Larson R. B., Tinsley B. M., & Caldwell C. M., 1980, 
ApJ, 237, 692 

Leitherer C, Schaerer D., Goldader J. D., et al., 1999, 
ApJS, 123, 3 

Leonardi A. J. & Rose J. A., 1996, AJ, 111, 182 

Leroy A. K., Walter F., Brinks E., et al., 2008, AJ, 
136, 2782 



41 



Lewis L, Balogh M., De Propris R., et al. 2002, MN- 
RAS, 334, 673 

Lucy L. B., 1977, AJ, 82, 1013 

MacArthur L. A., Gonzalez J. J., & Courteau S., 

2009, MNRAS, 395, 28 

Machacek M., Jones C, Forman W. R., Nulsen P., 
2006, ApJ, 644, 155 

Marcolini A., Brighenti F., & D'Ercole A., 2003, MN- 
RAS, 345, 1329 

McCarthy I. G., Frenk C. S., Font A. S., et al., 2008, 
MNRAS, 383, 593 

Mercurio A., Merluzzi P., Haines C. P., et al. 2006, 
MNRAS, 368, 109 

Merluzzi P., Mercurio A., Haines C. P., Smith R. J., 
Busarello G., Lucey J. R., 2010, MNRAS, 402, 753, 
Paper I 

Miller N. A., 2005, AJ, 130, 2541 

Mo H. J., Mao S., & White S. D. M., 1998, MNRAS, 
295, 319 

Monreal-Ibero A., Arribe S., Colina L., et al., 2010, 
A&A, 517, 28 

Moore B., Katz N., Lake G., Dressier A., Oemler A., 
1996, Nature, 379, 613 

Moster B. P., Maccio A. V., Somerville R. S. Naab 
T., Cox T. J., 2011, MNRAS, 415, 3750 

Moustakas J., Kennicutt R. C. Jr., & Tremonti C. A., 

2010, ApJS, 190, 233 

Murakami H., Baba H., Barthel P., et al., 2007, 
PASJ, 59, S369 

Nulsen P. E. J., 1982, MNRAS, 198, 1007 

Oh S., Oh K., & Yi S.K., 2012, ApJS, 198, 4 

Oosterloo T. & van Gorkom J., 2005, A&A, 437L, 
190 

Ocvirk P., Pichon C, Lancon A., Thiebaut E., 2006, 
MNRAS, 365, 46 



Peng C. Y., Ho L. C, Impey C. D., Rix H-W., 2010, 
ApJ, 139, 2097 

Pracy M. B., Kanschner H., Couch W. J., et al. 2009, 
MNRAS, 396, 1349 

Pracy M. B., Owers M. S., Couch W. J., et al., 2012, 
MNRAS, 420, 2232 

Proust D., Quintana H., Carrasco E. R., et al., 2006, 
A&A 133, 447 

Rich J. A., Dopita M. A., Kewley L. J., Rupke 
D. S. N., 2010, ApJ, 721, 517 

Rich J. A., Kewley L. J., & Dopita M. A., 2011, Apj, 
in press. 

Rickc G. H., Alonso-Herrero A., Weiner B. J., et al., 
2009, ApJ, 692, 556 

Roediger E. & Hensler G., 2005, A&A, 433, 875 

Roediger E. & Briiggen M., 2006, MNRAS, 369, 567 

Roediger E. & Briiggen M., 2008, MNRAS, 388, 465 

Rose J. A., 1985, AJ, 90, 1927 

Rossetti M., Ghizzardi S., Molendi S., Finoguenov 
A., 2007, A&A, 463, 839 

Sanchez S. F., Rosales-Ortega F. F., Kennicutt R. C, 
et al., 2011, MNRAS, 410, 313 

Sanchez-Blazquez P., Peletier R. F., Jimenez- Vicente 
J., et al., 2006, MNRAS, 371, 703 

Sanderson A. J. R. & Ponman T. J., 2010, MNRAS, 
402, 65 

Schlegel D. J., Finkbeiner D. P., & Davis M., 1998, 
ApJ, 500, 525 

Smith R. J., Hudson M. J., Nelan J. E., 2004, AJ. 
128, 1558 

Smith R. J., Lucey J. R., & Hudson M. J., 2007, 
MNRAS, 381, 1035 

Smith R. J., Lucey J. R., Hammer D., 2010, MNRAS, 
408, 1417 



Springel V. & Hernquist L., 2003, MNRAS, 339, 289 

Springel V., Di Matteo T., & Hernquist L., 2005, 
MNRAS, 361, 776 

Springel V., 2005, MNRAS, 364, 1105 

Steinhauser D., Haider D., Kapferer W., Schindler 
S., 2012, A&A, 544, A54 

Storey P. J. & Zeippen C. J., 2000, MNRAS, 312, 
813 

Sun M. & Vikhlinin A., 2005, ApJ, 621, 718 

Sun M., Jones C, Forman W., et al., 2006, ApJ, 
637L, 915 

Sun M., Donahue M., Roediger E., et al., 2010, ApJ, 
78, 946 

Sutherland R. S. & Dopita M. A., 1993, ApJS, 88, 
253 

Toomre A. & Toomre J., 1972, ApJ, 178, 623 

Tonnesen S. & Bryan G. L., 2009, ApJ, 694, 789 

Tonnesen S. & Bryan G. L., 2010, ApJ, 709, 1218 

Tonnesen S. & Bryan G. L., 2012, MNRAS, 422, 1609 

Urban O., Werner N., Simionescu A., Allen 
S. W.,B6hringer H., 2011, MNRAS, 414, 2101 

Valotto C. & Giovanelli R., 2004, AJ, 128, 115 

van Dokkum P. G., 2001, PASP, 113, 1420 

Vazdekis A., Sanchez-Blazquez P., Falcon-Barroso J., 
et al., 2010, MNRAS, 404, 1639 

Veilleux S. & Osterbrock D. E. 1987, ApJs, 63, 295 

Veilleux S., Cecil G., & Bland-Hawthorn J., 2005, 
ARA&A, 43, 769 

Veilleux S., Weiner B.J., Rupke D.S.N., et al., 2010, 
AJ, 139, 145 

Vogt N. P., Haynes M. P., Giovanelli R., Herter T., 
2004, AJ, 127, 3300 



Vollmcr B., Marcellin M., Amram P., et al., 2000, 
A&A, 364, 532 

Vollmer B., Braine J., Combes F., Sofue Y., 2005, 
A&A, 441, 473 

Vollmer B., Soida M., Chung A., et al., 2008, A&A, 
483, 89 

Vollmer B., Soida M., Chung A., et al., 2009, A&A, 
496, 669 

Vollmer B., Soida M., Chung A., et al., 2010, A&A, 
512, 36 

Wang Q. D., Owen F., & Ledlow, 2004, ApJ, 611, 
821 

Williams B. F., Dalcanton J. J.; Johnson L. C, et 
al., 2011, ApJ, 734, L22 

Yagi M., Komiyama Y., Yoshida M. et al., 2007, ApJ, 
660, 1209 

Yagi M., Yoshida M., Komiyama Y., et al., 2010, AJ, 
140, 1814 

Yoshida M., Yagi M., Okamura S., et al. 2002, ApJ, 
567, 118 



A Exploiting the simulations 

To compare the results of simulations with the data, 
the simulated data must be first correctly oriented 
with respect to the observer and then processed to 
mimic the characteristics of the observing instrument. 
These two steps are described in the next subsections. 

A.l Observing the simulations 

The output of the hydro dynamical simulations are 
tables with one row for each particle containing the 
three components of position and velocity. In our 
case the models are oriented with the spatial veloc- 
ity of the galaxy v ga \ along the (positive) Z-axis and 
the Y-axis on the galaxy disc. We first apply a ro- 
tation around the Y-axis to move the galaxy disc in 
the (X,Y) plane, so that \^ ga i hes in the positive XZ 
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semi-plane. Fig. [21] represents a coordinate system in 
which the disc of the galaxy lies in the (X,Y) plane. 
The direction to the observer is marked by the dashed 
lines ending with a small circle. Two viewing sit- 
uations, differing by a rotation of 180° around the 
X-axis, are in principle possible. In one case (case 
A) the wind points towards Z>0 (red arrows in the 
figure), while in the other (case B) the wind points 
towards Z<0 (cyan arrow). 

To project the model into the plane of the sky, we 
require that the spatial velocity of the galaxy projects 
into the observed line-of-sight velocity. The angle (ft 
between ~\? ga i and the line of sight is given by cos (ft = 
(Vlos/Vgai)- This is equivalent to require that wind 
makes an angle (ft with the line of sight in the direction 
of the observer. 

The top panels of Fig. [21] show the case where 
P =25° andV wind =2400kms- 1 , or (ft ~70°. The cone 
with /3=const. (hereafter '/3-cone'; blue) intersects 
the cone </>=const. (hereafter '0-cone'; black/gray) 
in two points in the Z>0 semi-sphere. The two black 
lines are the projection of the semi-planes passing 
through the Z-axis and the intersection points of the 
cones. They are at angles 6\ and 62 (—32°, 212° in 
this example) with the X-axis. Rotating the galaxy 
by one of these two angles, moves the wind direc- 
tion at the right angle with respect to the line of 
sight in which it projects into Vi os . In this case there 
are therefore two possible directions from which we 
can "observe" our model. In both views, the wind 
is directed as the red arrow in the upper-right panel, 
which corresponds to case A. Since the wind cone 
does not intersect the (ft-cone in the Z<0 hemisphere, 
no other viewing angles are possible. 

In the bottom panels of the same figure (corre- 
sponding to the case /3 =45° and (ft ~70°) the wind 
cone intersects the c/>-cone in both hemispheres. There 
are intersections in the Z>0 hemisphere at 61,62 (=- 
24.3°, 204.3°; black lines), and for Z<0 the intersec- 
tions take place at 63,64 (=-54.4°, 234.4° gray lines). 
Rotating the model by 6\ or 62, allows us, as before, 
to see the model from the correct point of view. But 
in this case the model can also be rotated around the 
Z-axis by #3 or 64. 

The second case corresponds to the cyan wind 
vector in the bottom-right panel and represents an 




Figure 21: Unit sphere (outline dotted), centred on 
the galaxy centre, with the X and Y axes on the galac- 
tic plane, and the Z-axis along the symmetry axis of 
the disc. The direction to the observer is marked by 
the small black filled circle in the semi-plane Y<0, 
Z>0, and is inclined by i=82° with respect to the 
Z-axis. In this way, the X-axis corresponds to the 
apparent major axis of the galaxy. The upper pan- 
els represent the case /3=25° and V lu j„d=2400kms~ 1 , 
while the lower panels represent the case /3=45° with 
the same Y w i n d- This wind speed corresponds to 
0=69.64°. The left panels offer a view of the unit 
sphere from the Z-axis. The blue circles are the in- 
tersection with the sphere of the cones /3=const. The 
black and gray ellipses are the intersection with the 
unit sphere of the semi-cone </>=const. directed to- 
wards the observer. The gray section is the half of 
the ellipse laying at Z<0, i.e. at the opposite side 
of the galaxy disc in this view. The right panels 
show the same configurations projected on the YZ 
plane. The dashed line is the line of sight. The dot- 
dashed ellipse in the upper-left panel is the intersec- 
tion of the unit sphere with the cone corresponding to 
V ru j n( i=1100 kms^ 1 , an example of a forbidden con- 
figuration. The arrows represent projections of the 
ICM wind (see text). 
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example of case B. To reach this configuration, the example is illustrated by the dot-dashed ellipse in 
model should be first rotated by 180° around the X- the top-left panel of Fig. EH which is the intersection 
axis. with the unit sphere of a 0-cone corresponding to 

The coordinate transformations needed in the above Y w i n d= 1100 km s _1 . This cone has no intersection 



process are the following. Let (3 be the angle of the 
wind with respect to the symmetry axis of the galaxy 
and i the inclination angle of the galaxy (plane of 
the sky to galaxy disc). First we move the symmetry 
axis of the galaxy into the Z-axis with a clock-wise 
rotation of (3 around the Y-axis. In this way, the 
wind lies on the (X,Z) plane. The second step is to 
move the wind vector in a direction in which the wind 
speed projects into \^ ga i- It can be shown that his is 
achieved by a rotation around the Z-axis by an angle 
9 given by 



cos 9 



cos 4> + k cos f3 cos 



sin (3 sin i 

where ft = — 1 for case A and k = 1 for case B. 
For each k this equation has two solutions: 9\ = 
arcos(cos 9) and 9% = tt — 9±. As a result, when 
only case A is allowed, there are two possible viewing 
angles (or 'projections'), while when both cases A 
and B are allowed, the possible projections are four. 

Other transformations are possible that do not 
change the physics of the models, but given that we 
fixed the position of the line of sight in the YZ plane, 
the only remaining meaningful transformation is a re- 
flection with respect to that plane. This corresponds 
to an inversion of the X-axis and to the change of sign 
of the corresponding velocity component of the model 
particles. This transformation makes the projection 
of the models on the sky to reflect with respect to 
the apparent major axis, and may become useful to 
match the models to the observed orientation of the 
galaxy, when this is not achievable with a simple ro- 
tation. In the present case this transformation was 
not needed. 

Cases A and B represent two quite different ob- 
serving conditions. In case A, the wind acts on the 
disc side opposite to the observer, in such a way that 
the outflowing gas overlays the disc. In case B the 
opposite is true, and part of the outflowing gas is 
hidden by the disc. 

Finally, we notice that not all combinations of 
parameters lead to an observable configuration. An 



with the /3-cone corresponding to (3=25°. This means 
that there is no possible orientation such that this 
configuration would lead to the observed line-of-sight 
velocity of our galaxy. 

A. 2 Simulating the observations 

The end product of the previous process is a matrix 
(xi,yi,V;; i=l,...,N partides ), where (x,y) are coordi- 
nates in the plane of the sky (x is along the apparent 
major axis), y is normal to x, V is the radial veloc- 
ity and ^particles is the number of particles in the 
simulation. 

To reproduce the effect of the seeing we substi- 
tuted each particle by a Gaussian distribution with 
the FWHM equal to the seeing. In each surface el- 
ement of the plane of the sky (see below), we com- 
puted the surface density of particles £ as the sum 
of the contributions from the Gaussians and the ra- 
dial velocity as the weighted average of the individual 
velocities. 

To speed-up the computation, we truncated the 
contribution of each particle at Wxacauss- This com- 
putation was performed on a grid of points spaced by 
0.1// on the sky, which we verified to offer a reason- 
able compromise between accuracy and computation 
speed. The points in the plane of the sky were even- 
tually sampled in 1// x 1// bins, matching the spatial 
binning of the data, and rotated to the position angle 
of the major axis of our galaxy. 
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